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We discuss the topological phenomena in the QCD-like theories with a variable number of fun- 
damental fermions Nf, focusing on the temperatures at or above the critical value Tc of chiral 
symmetry restoration. The nonzero average of the Polyakov line, or holonomy, splits instantons into 
(anti) self-dual dyons, and we study both the bosonic and fermionic interactions between them. The 
high temperature phase is a dilute gas of "molecules" made of 2Nc dyons, neutral in topological, 
electric, and magnetic charges. At intermediate temperatures the diluteness of the "molecular gas" 
reaches some critical level at which chiral symmetry gets restored: we comment on why it is different 
for the fundamental and adjoint fermions. At high density the ensemble is a strongly coupled liquid 
with crystal- like short range order: we speculate about its possible structure at small and large Nf. 
We finally show that certain lattice observations are in agreement with the proposed model, and 
suggest a number of further lattice tests. 



I. INTRODUCTION 

A. Qualitative overview 

Here we outline the qualitative findings that 
emerged during the course of this study, and 
provide the answers to some of the questions 
that followed from it. A history of the works 
and ideas that lead us to these answers will be 
provided in the next subsection. 

It is perhaps necessary to explain what we 
mean by the "dyonic vacuum." Dyons in gen- 
eral are objects possessing both electric and 
magnetic charges. However, the term is used 
in two very different contexts. Historically, the 
first one is the "particle dyons" (Julia-Zee), 
the excitations of 't Hooft-Polyakov monopoles, 
well known in the Georgi-Glashow model and 
many supersymmetric theories: we will not dis- 
cuss those in this work. The second type is 
the "self-dual dyons," which appear as con- 
stituents of the instantons. Rather than be- 
ing particlelike excitations of the vacuum, as 



the monopoles and the first type of dyons are, 
they are part of the vacuum itself, describing 
a certain topologically nontrivial configurations 
of the gauge fields. Not being particles, they do 
not have momenta or kinetic energies: they ap- 
pear in the QCD partition function integrated 
over their collective variables. While instantons 
had among such variables color orientations, 
the dyons have their positions and Abelian (di- 
agonal) color charges. In the case of the SU(2) 
color group, there are two dyons per instan- 
ton, commonly called the M-type and the L 
(or twisted) type, (see Table [l|. For a review 
of instanton dyons see [l], p and references 
therein. 

Here are some physical questions we will dis- 
cuss: 

1. What are the interactions between dyons, 
especially between the self-dual and anti- 
self-dual ones? 

2. How do fermions contribute to the inter- 
action between dyons? 
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3. What is the quahtative picture of the 
dyon statistical ensemble, as a function of 
three key parameters, the temperature T 
and the number of fundamental Nf or ad- 
joint Na quarks in the theory? In partic- 
ular, why does the chiral transition moves 
to a stronger coupling? 

4. In the high-T limit, gauge field topology 
was described as a dilute gas of instanton- 
anti- instant on molecules [3 . How are 
these objects modified for the case of 
the nonzero holonomy, in the language of 
dyons? 

5. What is the Dirac eigenvalue spectra for 
different dyonic ensembles? At which 
condition chiral symmetry breaking takes 
place? 

6. Can one explain the dependence of the 
chiral phase transition on Nf and/or Na^. 

7. Can one evaluate the "gaps" in the Dirac 
eigenspectra which are developed at T > 
Tc using the dyonic ensemble? 

8. At T > Tc, using the quenched ensemble 
of gauge fields, it has been found on the 
lattice that the chiral properties crucially 
depend on the particular periodicity con- 
ditions for the fermions. In particular, 
the periodic ones do not show a chiral 
restoration transition, unlike the (physi- 
cal) anti-periodic fermions. How can one 
understand these observations? 

9. Why does the chiral transition strongly 
depend on the color representation of the 
fermions, such as the fundamental or ad- 
joint ones? 

The reader perhaps noticed that this list of 
questions includes neither a discussion of the 
holonomy potential nor other questions related 
to confinement (such as e.g. in [36 .) We think 
that any assessment of the back reaction of the 
dyons on the holonomy can only be done after 
a more quantitative theory of their ensemble 
emerges. The purpose of this paper is to take 
a step toward developing such a theory. Thus 
here we will consider the holonomy (P(T)) as 
given, e.g. by the lattice data. 



Now we outline the picture. It is convenient 
to discuss it by defining three regimes, from 
high to low T. We will call them: 

(i) High T ^ Tc case, the regime of a dilute 
molecular gas 

(ii) Intermediate regime, T = (1..2)Tc, inter- 
acting molecular gas 

(iii) Dense regime, T < Tc^ dyons form a 
strongly coupled plasma, in their liquid 
phase 

Here are further explanatory comments on 
each of those: 

(i) High temperature implies weak coupling 
and thus the semiclassical treatment of instan- 
tons/dyons is applicable. Since these objects 
have nonzero electric fields, subject to the per- 
turbative Debye screening, their density at high 
T is strongly suppressed. In a resulting di- 
lute regime, the ensemble forms a "molecu- 
lar gas" of objects that have zero topologi- 
cal, electric, and magnetic charges. The av- 
erage Polyakov line in this regime is close to 
1, i.e. the "Higgs Vacuum Expectation Value 
(VEV)" ~ 0, so all M-type dyons are light, 
while "heavy" L, L dyons have nearly all the 
action of instantons/antiinstantons. Instanton- 
antiinstanton molecules were described in 
and subsequent works: and in the high-T limit 
we expect to be close to those results. 




FIG. 1: (Color online) The schematic picture of 
the dyonic molecule, for 2 colors and large Nf. 



Fermionic zero modes of the instantons are 
shared by their constituents in a way that de- 
pends on the type of fermion of the theory. 
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Physical antiperiodic fermions have zero mode 
of the (twisted) L, L dyons. As the number 
of fundamental fermions Nf in the theory in- 
creases, they bind them into tight LL "clus- 
ters," which play a role of the nucleus of these 
molecules. Consequently, the chiral symmetry 
is unbroken and the lowest Dirac eigenstates "at 
the gap" correspond to independent LL clusters 
(see Fig. [l]). 

The standard Abelian electric charges of both 
L and L are equal to -1, so the clusters has 
the charge -2. (The molecule thus looks like 
anti-He, with M, M as "positrons" around it). 
A particular sign of a charge does not violate 
C parity, of course, because the Abelian fields 
are color-projected to the color direction of the 
holonomy field ~ Tr (^^^^^4), or more precisely 
~ Tr F^yL^ where L is the Polyakov loop. Since 
both non-Abelian fields in it are C odd, the 
product is C even. These signs are just matter 
of the definitions used in the field. 

Lattice practitioners sometimes use the so 
called "valence" or non-dynamical quarks (not 
included in the partition function) as a tool for 
the analysis of the gauge configurations. Those 
may have arbitrary properties and periodicity 
conditions. The "valence quarks" periodic over 
the Matsubara time have completely different 
zero modes and interact with the lighter M- 
type dyons. Those are also more weakly corre- 
lated than the L-type ones. The difference in 
their Dirac eigenspectra will be important tool 
in testing the structure of the dyonic vacuum. 

(ii) As T is lowered, the effective electric cou- 
pling as{T) = ^^(T)/47r grows and eventually 
becomes large as = 0(1). A quite specific point 
introduced in [5 is the so-called "E/M equilib- 
rium point," at which as{T) = 1. Because of 
the Dirac condition for electric and magnetic 
couplings 

Q^sQ^magnetic = integer (1) 

at this point, with integer being 1, the magnetic 
alpha is also 1. It has been argued and con- 
firmed on the lattice that at this point magnetic 
excitations - monopoles - become about as nu- 
merous as the usual electric excitations, quarks 
and gluons. In gluodynamics (no fermions, 
Nf = 0) this happens at T ^ 1.4 Tc: how this 
depends on the presence of fermions remains to 
be studied. 



Below Tc the confinement forces the density 
of the electric objects (quarks and gluons) to 
be zero, while the magnetic (monopoles) retain 
the finite density. Only at T ^ does it disap- 
pear as well, with only the magnetic condensate 
remaining. Consequently, the electric Debye 
mass disappear at T < Tc, while the magnetic 
screening mass remains finite. This implies that 
the electric screening of instantons at high T is 
substituted by a magnetic screening at T < Tc. 
As argued in [38 , the latter creates a factor in 
the density of instantons ~ exp(— x const.), 
where p is the instanton radius and the con- 
stant, due to magnetic screening, has a nonzero 
value even at T ^ and is related to Base- 
condensed magnetic monopoles. In certain dual 
models this constant was further related to the 
QCD string tension const = 27ra [38]. The ex- 
pression describes well the lattice data on in- 
stanton size distribution and also explains why 
in the QCD instanton ensemble remains rela- 
tively dilute even at T = 0. It would be very 
interesting to see if any of that remains to be 
true at large Nf. 

Near Tc magnetic screening should be in- 
duced mostly by the scattering of non-Bose- 
condensed monopoles. To our knowledge no 
study of the effect has ever been done, and we 
also defer it to future studies. 

The interrelation between the "particle- 
monopoles" (inducing confinement by their 
BEC) and the instanton-dyons (inducing chi- 
ral symmetry breaking as we discuss below) 
is of course an intriguing open problem. On 
the level of gauge configurations themselves or 
their zero modes one finds no direct relation 
between them. However, at the level of the ef- 
fective 't Hooft Lagrangian an intriguing rela- 
tion has been found by Poppitz and Unsal ^] 
in N=2 compactified Super Yang Mills (SYM) 
case. It appears at the level of the partition 
functions, one being the sum of the particle- 
dyon excitations and another the sum over the 
periodic instanton-dyon semiclassical solutions. 
Such relation clearly deserves further study and 
generalizations. 

While any perturbative expressions/intuition 
is obviously not to be trusted in the regime 
with as{T) ^ 1, the lattice simulations treat 
this region consistently. Furthermore, as we 
will detail below, in QCD-like theories with 
many fermions the plasma phase extends to 
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even stronger coupling of ~ 3 or so. Perhaps 
the dual-magnetic-formulation of such theories 
can be used in this case, as the magnetic cou- 
pling is getting weak amagnetic ^ 1/3. 

Our paper, as many others, rely on the ro- 
bustness of the topological effects under defor- 
mations, even if the amplitude of those is not 
small. Furthermore, topology is related (by 
index theorems) to fermionic zero modes. A 
"collectivized" set of such fermionic states con- 
tributes significantly to the quark condensate, 
pions and strongly influence the structure of the 
lowest hadronic states and correlation functions 
[7 . In contrast to the papers discussed in that 
review, we now approach this problem "from 
above," starting from the hot symmetric phase 
at T > Tc and follow the evolution of the topo- 
logical quark states, from localized to delocal- 
ized ones as the transition temperature is ap- 
proached. 

Lattice data tell us that in the tempera- 
ture interval (2..1)Tc the average Polyakov line 
(P(T)) changes from ^ 1 to ~ 0. The holon- 
omy changes from i^ = 0toi^ = l/2, at which 
point the masses (i.e. actions) of the L,M-type 
dyons become comparable. We also know that 
at the latter point confinement phenomenon 
takes place. 

(iii) We can only qualitatively discuss the 
dense regime near and below Tc, as the inter- 
action between the dyons gets very strong. We 
try to approach the problem from the perspec- 
tive of the strongly coupled classical plasma. 

For large Nf the basic objects include the 
LL clusters which, we will argue, strongly repel 
each other. Therefore, the optimal correlations 
in such a medium would be similar to other sys- 
tems which experience strong repulsive forces, 
such as closely packed liquids. While the global 
order is absent, locally those are strongly corre- 
lated, with the type of correlations being simi- 
lar to those in certain best-packed crystals. 

For zero Nf the dominant forces are 
Coulomb-like and corrections to them are in a 
form of the determinant proposed by Diakonov 
and collaborators, as well as the screening ones. 
If so, we suggest dyonic crystals resembling salt: 
cubic with alternating L, M dyons. 

(Needless to say, we do not think that the 
solid phase is actually reached; it is well known 
that strongly correlated liquid have short-range 
correlations that are the same as their fully 



ordered, crystalline form. While the "dyonic 
crystals" discussed provide examples of config- 
urations in which the interaction is minimized, 
thermal fluctuations do kill the long range or- 
der, making it a liquid. Perhaps it is worth 
mentioning that the main parameter in the 
Boltzmann exponents, the mean ratio of the in- 
teraction potential per particle to T, also called 
r, needs to be F > Fc ~ O(IOO) for solidifica- 
tion. In the dyon problem discussed this is not 
so, as F - 0(10) < Fc.) 

Let us now focus on the main observable 
to be discussed, the Dirac eigenvalue spectrum 
and possible chiral symmetry breaking. With 
the increasing number of fundamental fermions 
Nf in the theory, they induce stronger corre- 
lations and reduce the size of the LL "clus- 
ters." If one wants to follow the lines of con- 
stant quark condensate, e.g. the chiral restora- 
tion line, one has to increase the density of the 
clusters accordingly. This can only be achieved 
by a shift to stronger coupling. As the dyon 
masses and interactions are ~ 1/^^, they be- 
come lighter and less interacting. (Needless 
to say, their fermionic zero modes and related 
interaction must still be there, for topological 
reasons: they do not depend on the coupling.) 
The (dimensionless) density of dyons continue 
to grow to the situation in which the inter- 
dyon distances become comparable to the LL 
molecule size. 

The adjoint fermions are very different from 
the fundamental ones. Some of them remain 
"massless" (in the sense of the "holonomy 
mass") after adjoint Higgsing, and this dras- 
tically changes the dependence of the "hopping 
amplitudes" on the distance, from exponential 
to powerlike. The chiral symmetry is unbroken 
above such T when not only heavy L but also 
light M dyons have zeromodes. This puts the 
chiral phase transition at much weaker coupling 
(higher T) 

B. Prom instantons to dyons 

The discovery of the instanton solution [6] 
has created a great deal of literature, including 
electroweak physics of baryon charge noncon- 
servation as well as the famous exact results for 
various super symmetric theories. Obviously we 
cannot review this amount of literature here. 
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In the context of the QCD-hke theories, the 
predecessor of this paper is the so-cahed instan- 
ton hquid model, for a review see [7 . Its main 
point was to account for the so-cahed 't Hooft 
interactions to arbitrary order, by including 
the fermionic determinant in certain approxi- 
mation in numerically simulated statistical en- 
sembles. The calculated point-to-point corre- 
lation functions have reproduced many lattice 
results related to chiral SU{Nf) and U{1) sym- 
metries. Chiral restoration can be viewed as 
the disappearance of the nontrivial solution to 
the so called gap equation. Alternatively, it 
was explained [3 as a consequence of a struc- 
tural phase transition in the instanton ensem- 
ble, from a random plasma at low T into a gas of 
strongly correlated // instanton-anti-instanton 
pairs. The pairing mechanism is due to the 
fermion exchange, thus it gets stronger as Nf 
grows. 

Let us recall its basic ideas which will be 
used below. In the basis spanned by the zero 
modes of individual instantons/anti-instantons, 
one can write the Dirac operator as 



TiA 
Tai 



(2) 



where we have introduced the overlap sub- 
matrix TjA 



TiA 



(fx^l)l j{x- zi)ipilJo,A{x - zaI^) 



where /, A are indices which run over all in- 
stantons and antiinstantons in the configura- 
tion. Here, ipoj is the fermionic zero mode. 
The individual matrix elements have the mean- 
ing of a hopping amplitude for a quark from one 
pseudoparticle to another, and the determinant 
of this matrix is nothing else but the sum over 
the loop diagrams in which quarks visit each in- 
stanton once. Note that two ?/^s have opposite 
chirality, so if i=instanton then j=antiinstanton 
or v.v. The fermionic determinant is approxi- 
mated by \det{Tij)\'^ . The low-T ensemble is a 
dense liquid that breaks chiral symmetry, but 
at high T (small size in r direction) it breaks 
into "// molecules" and chiral symmetry gets 
restored. The actual calculation was a simula- 
tion of the ensemble with the weight containing 
\det{Tij)\'^ , which was then used for the evalu- 
ation of the Dirac spectra and hadronic corre- 
lation functions. At high T the approximate 



factorization of the Dirac matrix into indepen- 
dent 2x2 boxes (for separate clusters) explains 
the deformation of the Dirac eigenvalue spectra 
and disappearance of near-zero eigenvalues and 
the existence and the magnitude of the spectral 
gap G. 

We will now extend these ideas to the case 
of the nonzero holonomy, the gauge- invariant 
closed loop integral over the = r cir- 
cle f^^^ drA^. Its exponent, the so called 
Polyakov line, averaged over the statistical en- 
semble of fields, has a nonzero value 



Tr exp 



drA^ ])^0 (4) 



This calls for classical solutions that do not ap- 
proach zero fields at spatial infinity but rather 
some constant value v ofthe A^ (in SU(2)). We 
will also use dimensionless notations 



27rT' 



(5) 



Explicit solutions of such type [H [TO] demon- 
strate that an instanton gets split into the Nc 
constituent dyons. The names and quantum 
numbers (for the simplest SU(2) gauge group 
we will discuss in this work) cover all four pos- 
sibilities for the electric and magnetic charges, 
see Table [l[ For SU {Nc) in general there are 
Ml, M2...Mnc-i static dyons with all diagonal 
charges and one "twisted" L-dyon. 



name 


E 


M 


mass 


M 


+ 


+ 


V 


M 


+ 




V 


L 






27vT -V 


L 




+ 


27iT-v 



TABLE I: The charges and the mass (in units of 
Sn^/g^T) for 4 SU(2) dyons. 

Let us indicate here the qualitative differ- 
ence that the nonzero holonomy brings into 
this problem. The fundamental fermions in 
the "Higgsed" vacuum with a VEV of A4 are 
"massive," [50 , with masses (in SU(2)) m/ = 
:^gv/2. Therefore, the zero modes at large dis- 
tances r ^ 00 decrease exponentially with the 
distance, unlike the power behavior typical for 
the zero holonomy case. These rapidly decreas- 
ing fermionic amplitudes are of course further 
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enhanced by the number of fermionic zero 
modes 

e-^-detT-e-^-^^^^ (6) 

which creates strong hnear confining potential 
for the corresponding dyons and thus produces 
smah-size "clusters" of the size 

(r) ~ (N^nif)-' (7) 

The number of the modes is dependent on 
the fermion's color charge and the number of 
its copies. For the usual fundamental quarks 
= 2Nf^ as there is a zero mode for a quark 
and for an antiquark. 

For the adjoint fermion = 2NcNa. Fur- 
thermore, adjoint fermions that are diagonal 
with respect to the Polyakov line VEV remain 
"massless." The interaction between the dyons 
due to an exchange of the adjoint fermions has 
been discussed by Unsal [TTj . The center of his 
proposal is a "bion," a cluster of LM dyons 
with the magnetic charge 2. Such bions can in- 
duce Polyakov-like confinement in the spatially 
compactified QCD pjj. 

C. Overview of lattice data on chiral 
symmetry restoration and deconfinement 

1. The critical lines versus the number of 
fermions Nf , Na 

Let us start by reviewing some recent lattice 
results. Our version of the phase diagram uses 
the "critical lattice coupling" 

Ara = ^ (8) 

as a function of Nf or Na. The "bare" coupling 
values in lattice works are defined at the lattice 
UV scale a. In order to make it lattice inde- 
pendent, we have evolved the scale from a by a 
factor of Nt (the number of points in the time 
direction) to the physical scale Nta = l/Tc us- 
ing the two-loop beta function. The near over- 
lapping points in Fig. |2|a) are from different 
Nt simulations: their spread is a measure of the 
inaccuracy of the two-loop beta function used. 

The open diamonds in Fig|2|a) correspond 
to lattice data from Ref. [13^ : they show 



the critical line of the chiral restoration (thin 
solid line). Above the line one finds symmetric 
[quark-gluon plasma (QGP)] phase, while be- 
low is the chirally broken (and confining) one. 
Moving downward on this figure means increas- 
ing the gauge coupling or decreasing the tem- 
perature. Thus, dit Nf ~ 10 one may reach 
"the most strongly coupled QGP," which is by 
itself a very interesting phenomenon. Why is it 
happening? 

(There are many other simulations reported 
in the lattice literature, of course: we decided 
not to put those on this plot because the rather 
different actions used produce rather random 
spreading of the couplings, confusing the pic- 
ture.) 

The situation at Nf = 12 is special. Ref- 
erence [H] argued that the chiral symmetry 
remains unbroken. A more recent paper 
have studied the region of even stronger cou- 
pling ~ 1, and perhaps clarified the pic- 
ture. Two distinct phase transitions are ob- 
served, with the chiral condensate disappear- 
ing at stronger coupling (the closed diamond 
in Figj2|a)) than the confinement (the closed 
box). Thus, a novel intermediate phase in be- 
tween is confining but chirally symmetric. An 
understanding of its precise nature remains a 
challenge, although such examples are known 
in super symmetric theories. 

The same authors [15 also concluded that 
for Nf = 12 those transitions are also separated 
by a bulk transition from a weaker coupling do- 
main, in which there seem to exist a conformal 
(infrared fixed point) behavior. If so, the two 
phase transitions must be below (on the other 
side of) the true trajectory of the fixed point 
schematically shown by the (red) dash-dotted 
line. 

[This line may deviate from the thick solid 
(red) line, showing a fixed point line using the 
two-loop beta function, because of the very 
strong coupling involved. The issue of where 
exactly the conformal window starts remains 
unresolved. We tentatively put the vertical 
dashed line separating it at A^/ = 11 in the fig- 
ure. This is not a prediction but just a guess, 
and should not be used in any way.] 

Changing the fundamental quarks into either 

(i) adjoint (triplet in SU(2), octet in SU(3)) or 

(ii) symmetric rang-2-tensor (triplet in SU(2), 
sextet in SU(3)) shifts the Tc of the chiral tran- 
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FIG. 2: (Color online) The critical lines for chiral restoration (solid line and diamonds) and deconfinement 
(dotted line and boxes) of the Nc = 3 gauge theory. We plot the critical lattice coupling /3c{Tc) — Q/g1{Tc) 
versus the number of fundamental quarks Nf in (a) or a number of adjoint quarks Na in (b). Both paths of 
the infrared fixed point, calculated in the two- loop approximation, are shown by the thick (red) lines. The 
vertical dashed lines separate the "conformal window domain": its location is a guess. In Fig. (a) we also 
show, by the dash-dotted (red) line, our guess for the actual location of the fixed point. For the meaning 
of the data points see the text. 



sition upwards. The theory with two adjoint 
fermions, Na = 2, has been studied in detail, 
see Fig. ^h) based on [16 , but (to our knowl- 
edge) not for other values, notably Na = 1. 
(Of course, introducing variable masses may 
allow to follow the lines continuously.) Two 
distinct transitions were found, but in the op- 
posite order to the previous case of the fun- 
damental fermions. So, between the solid and 
the dotted lines there exists the deconfined but 
chirally broken phase (a plasma of constituent 
quarks). Furthermore, while the difference be- 
tween these two points may appear small on 
this plot, the actual Tc scales of the two transi- 
tions are different by about factor 8. 



(It is also noteworthy that these points are 
close to the conformal window perturbatively, 
or perhaps already inside it nonperturbatively. 
We tentatively put its boundary - the vertical 
dashed line- at A/'a > 2: it is not a prediction. 
Its exact location needs to be found numerically 
in future works. ) 



2. The magnitude of the chiral splittings versus 
Nf. 

For Nf up to at least 8, the low-T theo- 
ries retain both confinement and chiral symme- 
try breaking, but the quantitative relation be- 
tween them changes. Let us characterize it by 
the relative splittings between the chiral part- 
ners, such as vector-axial p — Ai mesons or the 
nucleon and the lowest l/2~ A/"* resonance at 
T = 



TUN 



(9) 



(10) 



the Nf dependence of these ratios is interesting: 
in the interval Nf = 0..3 these chiral splittings 
are "large" near the experimental values (so to 
say, at Nf ^ 2.5) 



Aexp 



0.45 



xexp 



0.55 



(11) 



which are well reproduced on the lattice. Yet 
calculated for Nf = 4 ^/Tj and 8 \1S^ theories 
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one consistently finds about twice smaller val- 
ues, and at Nf = 12 these splittings were not 
observed at all [18^ . In view of the trend 
just discussed, as well as because of the two 
transition observed in [15] , we expect the chiral 
and deconfinement lines to separate, perhaps 
already at A^/ > 4 or so. 



3. Sensitivity to the fermionic periodicity 
conditions. 

Introducing an arbitrary phase in the peri- 
odicity condition of the "valence quarks", one 
can switch the fermionic zero mode between the 
dyons: this has been demonstrated using arti- 
ficial configurations for calorons, e.g. in |19] . 
There is significant literature covering efforts 
efforts to understand the difference in lattice 
gauge configurations below and above Tc. A pa- 
per presenting interesting results on the Dirac 
eigenvalue spectrum in the SU{3) quenched 
and unquenched ensembles is that by Bilgici et 
al [20 . Its brief summary: 

(i) at T > Tc the Dirac spectrum has a well- 
determined gap G (no eigenvalues inside A < 
|G|), growing approximately linearly 

G^{T-T,),T>T, (12) 

in the quenched case, till at least about 2Tc. 

(ii) if arbitrary (twisted) boundary conditions 
are used for (valence) fermions, by a phase 
(j) = 2tiz in a periodicity condition, they seem 
to be irrelevant below Tc but change the re- 
sults drastically above it. {i^il))^ or th density 
of eigenvalues at zero, seem to have a simple 
dependence on the angle 

KV^V^)! ^ci(T)+C2(T)cos(/) (13) 

with only one harmonics and positive coeffi- 
cients ci,C2. For the holonomy values shifted 
above Tc by ±27r/3 the phase of the cosine is 
shifted accordingly. 

(iii) As a result antiperiodic fermions cos0 = 
— 1 have a density touches touches zero at Tc 
and develops a gap, restoring chiral symmetry. 
The periodic fermions, with cos(/) = 1, never do 
so and their ('^/SV^) grow indefinitely above Tc. 



4. Chiral restoration in different Polyakov phase 
sectors. 

The previous issue is strongly related 
to lattice observations that subensembles of 
quenched configurations at T > Tc with dif- 
ferent Polyakov phases show different spectra 
of the Dirac eigenvalues and chiral parameters. 
In SU(3) there are two sets, one with (P) real 
and another with the phase e^*^^/^. For exam- 
ple. Fig 1 of [22 shows that these two different 
subensembles not only have different eigenval- 
ues but also drastically different "participation 
ratios" (a degree of mode localization on the 
lattice). This phenomenon gives us the oppor- 
tunity to study more than one set of holonomies 
in one simulation. 



5. The spectral gaps at T > Tc versus the 
fermion periodicity conditions. 

The gap opening is clearly observed for an- 
tiperiodic fermions but not for periodic ones, 
see [20 . An explanation based on dyon- 
antidyon classical correlation at high Tc was of- 
fered in [4]. Similar studies based on quenched 
configurations have been extended to adjoint 
fermions in Ref. [21]. Like for fundamental 
quarks, the antiperiodic adjoint fermions show 
a clear gap opening above the chiral transition, 
while the periodic ones do not. 



6. Identifying the topological objects via their 
fermionic zero modes. 

In configurations with say exact Q=l one 
can locate one exactly zero mode and see 
the location of the corresponding eigenvector 
is space- time. Using such lattice configura- 
tions Gattringer and Schaefer [23] have ob- 
served that while the eigenvector does indeed 
locate a single "topological lump," its position 
and quantum number depends on the bound- 
ary phase and jumps at certain values, resem- 
bling what happens with the different types of 
constituent dyons inside classical caloron solu- 
tions. Such techniques allow locating all kinds 
of dyons and potentially study their correla- 
tions / interactions. 
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At T ~ Tc, with massless fermions and re- 
stored chiral symmetry, all configurations with 
nonzero Q are absent, and the Dirac eigen- 
values get gapped. What are the states "at 
the gap," with the lowest Dirac eigenvalues at 
T > Tc^ As demonstrated in [24 , those have 
two topological lumps, confirming the picture 
of the paired instanton-antiinstantons [3 . As 
we will argue below, in the language of dyons, 
these molecules are more complicated, with 2Nc 
dyons of all kinds and certain Abelian charge 
distribution. Therefore, now one should use 
this method again, to look at Abelian-projected 
charges, clarifying their structure further. 

More recently, Bruckmann et. al. [25] looked 
at the fermionic states at the chiral gap at 2.6 
Tc for quenched SU(2) gauge simulations. They 
observe that the corresponding eigenstates are 
well localized and correspond to a strongly 
modified local value of the Polyakov line. They 
have shown that the number of such objects 
vastly exceeds the density of the isolated topo- 
logical charges deduced from topological sus- 
ceptibility, ruling out an ideal instanton gas as 
their source. They also commented, at the end 
of the section on topology, that their data "do 
not exclude" configurations in which the topo- 
logical charge cancels, like the instanton-anti- 
instanton molecules to be discussed. 



II. DYON INTERACTIONS 
A. Classical Interaction 

As is well known, classical interaction of the 
dyons inside one of the sectors - self-dual or 
anti-self-dual - are absent, as they are protected 
by the Bogomolny-Prasad-Sommerfeld (BPS) 
bound. Although it is clear how this works from 
the explicit solutions [9l[T0], we will still discuss 
it in the dilute limit, as our starting point. 

The question of Higgs topology and 
monopole interaction has been addressed 
in many referances (e.g. |45H47j ). Recall 
that the usual Higgs field has to go to zero 
at the monopole center because there is no 
preferred direction of color there. However 
with our "Higgsing" by the Polyakov loop, 
which is an element of 5/7(2), and its value 
at infinity [which defines an unbroken U{1) 



direction], the Polyakov loop can be viewed 
as L G SU{2)/U{1) = S'^ mapping. At points 
without definite color direction the Polyakov 
loop takes a value of L{xo) = ±1: thus, two 
types of dyons. Indeed, the effective "Higgs" 
field at the dyon's center it is going to zero for 
M type and to A4 = 7T(p ■ r for the "twisted" 
L dyon, where (/) is a unit color vector (see 
below). 

For a caloron -the L — M pair in SU(2) - 
one can see from quantum numbers that both 
the electric and magnetic charges are opposites, 
so they should both create attraction. Another 
long-distance force is Higgs mediated: because 
dyons of M and L types have different value 
of the A4 at their centers this turns out to be 
repulsive. Furthermore, it exactly cancels the 
attractive Coulomb forces, as required by BPS. 

Now we consider dyon-antidyon pairs, start- 
ing with MM (which do not have a temporal 
twist, i.e. dyons completely static in time). We 
take them at a large distance d » 1/v com- 
pared to the size of their cores. Inside some 
ball around the dyon (antidyon) of radius Tq 
such that 1/v « ro « d the field strength 
can be written as a small deviation from self- 
duality due to the other dyon, i.e. of the order 
1/d . Outside of these balls the "Higgs" field is 
given by {A^l = v — 1/ri — l/r2, where ri,2 are 
distances from the dyon and antidyon, so as to 
conform to the expected asymptotic formulas 
for dyon and antidyon (see e.g. [2 ). 

For a single dyon the Higgs field can be writ- 
ten as 

A4 = h{r)4> , (14) 

where (j) = A^/IA^]. Asymptotically h{r) = 
V — 1/r. An ansatz that properly describes two 
dyons would have to obey the condition that 
asymptotically 1^441 ^ v — l/ri — l/r2. However 
the color direction is a gauge choice that can be 
chosen arbitrarily at each point. Insisting that 
the Higgs points are in one color direction at 
some large sphere (gauge combing), one then 
has to introduce Dirac strings, as the gauge 
transformation cannot be made single valued. 
We do not specify this gauge choice, as we only 
deal with the action, which is gauge invariant. 

That being said, we expect that the influence 
of the other dyon will change the h function by 
introducing an additive Coulomb term near the 
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core of the first dyon, i.e. if ri << r2 we have 



F(ri,r2)«(%i)-l/r2) 



(15) 



where ri and r2 are distances from monopoles 
to the point of observation. An analogous rela- 
tion holds when r2 << ri. 

Now we determine the action of the system 
of two dyons, writing the action as an integral 
over three regions 

5 = ^ (^ TV {F^)d^x + {F^)d^x 

+ f Tr {F^)d^x\ (16) 

J outside / 

where 1 and 2 denote the regions around the 
dyon and antidyon respectively. Inside these re- 
gions we assume self-duality (anti-self-duality) 
up to some small correction of the order of 1/d 
inside the cores, i.e. DiA^ = T^^ijkF^^ -\- fi, 
where fi is the field strength deviation from 
self-duality induced by the other (anti)dyon 
and is of order f = o{l/d). 



/3 J (Tx di{A^B') + o(l/d^) = 
4.(.-l-i)+o(l/d^), 



(17) 



where the integration is over a ball of radius ro 
centered around the first dyon. 

We pause to comment on similar expression 
in the case when we have a single dyon. One 
can integrate on the surface at infinity, which 
would just simply yield AiivP, i.e. the usual 
mass of a dyon. However, it is more instructive 
to divide the region of integration into a small 
ball of radius ro and the rest. The small ball 
is a total derivative and yields the contribution 
47r('U — l/ro)/3) to the action. Then we can write 
the action as 

S single dyon = ^7r{v - l/ro)/3 (18) 
J outside ^ 

where we integrate over the volume outside 
the ball. However in this region the fields are 



Abelian and behave in an expected way. We 
can write Ei = ^^(j) and Bi = and the 

outside integral is 47r/ro. The sum of the re- 
gion inside and outside gives the expected re- 
sult Attv. 

Coming back to the case of two dyons, we 
include the region around the antidyon and get 



1. 



-Tr 



/ 



d^x = Stt 



1 

^0 



(19) 



showing how one dyon has been modified by 
the presence of the other one. Next we write 
the integral outside as the sum of electric and 
magnetic parts, i.e. 



E,. = 



J. _n_n 

^i 9 9 



(20) 
(21) 



which, upon integration, give the expected in- 
teraction 4:7r/d for the electric and —4:7r/d for 
the magnetic, and they cancel. Also there are 
two self-energy terms which are given by 47r/ro 
for electric and magnetic field separately, which 
cancel the 1/ro contribution to the inside-the- 
sphere integration. 

Notice that even though the electric and 
magnetic fields cancel outside the cores, the 
dyon-antidyon system still attracts, due to the 
modification of their mass by the presence of 
the other (anti)dyon. Thus there exists a long- 
distance classical Higgs-based attraction for the 
MM pair. 

One can equally well consider LL dyons with 
the twist, i.e. with core time dependence. The 
only modification is that then the contribution 
to the action of the core is given by Airv, where 
V = 27T — V. Also for a purely self-dual sector, 
the interaction of the L and M dyon is seen to 
cancel. This result is well known from the orig- 
inal works Ho] , but here we see that since the 
"Higgs" asymptotic looks like v — I/vm + 1/^l 
(see [2 ), the "Higgs" interaction is repulsive, 
which exactly cancels the attractive forces of 
the E and B field. Note also that due to this ef- 
fect of "dyon mass renormalization" we expect 
that if the L dyon has a fermionic zero mode 
(which as we will see in the next subsection, 
depends on this holonomy), it is renormalized 
by the presence of the other M dyon. This was 
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observed in the original papers by van Baal et. 
al. and it followed from the exact zero mode 
expression: we just identify its physical origin. 



B. Fermion- induced interactions 

Fermionic interactions between dyons are 
central for this paper. They are induced by 
the presence of fermionic localized modes facil- 
itated by the time-dependent twist in the gauge 
fields. The fermions introduce the fermionic de- 
terminant det^ factor in the partition func- 
tion. If both the dyon and anti-dyon are in 
isolation (at large distances), they have zero 
modes, which leads to a vanishing determinant: 
thus, such configurations are excluded from the 
ensemble. Obviously, at finite r the modes are 
nonzero, and therefore the dyon-antidyon pair 
is attracted due to the fermions. 

As was done for the instantons, one can look 
at the Dirac operator in the basis of the lo- 
calized zero modes of the individual L and L 
dyons. The matrix element of the p zero mode 
between two of those we denote as T/j, where 
the indices run through all dyon and antidyon 
zero modes. Since the Dirac operator in the 
chiral basis connects between the left and right 
fermions only, the diagonal elements are all 
zero, and only blocks Tjj and Tjj remain, where 
now / runs through the dyon zero modes only, 
and / through antidyon zero modes only. 

It is quite clear that in the case of a dyon- 
antidyon pair, since detT = — we have 
that Veff = — ln(T^i;(r^/;)). Since the matrix 
element is approximately Tj^i ~ g-Mr/2^ where 
M is the "holonomy mass" of the fermion, the 
resulting effective potential between dyons is 
linearly confining. 

While in this paper we will focus only on 
the zero modes of the fundamental quarks, we 
would like to mention some important works on 
the adjoint fermions, which naturally appear in 
the supersymmetric context. For periodic com- 
pactification the corresponding index theorem 
is discussed in [39 (see also citations therein). 
An extensive discussion of the zero modes for 
the periodic and antiperiodic adjoint fermions 
can be found in [40] • 



1. Fermionic zero mode for arbitrary periodicity 
condition 

As we mentioned in the Introduction, one 
should not confuse the "particle dyons" and 
"instanton" (or self-dual) dyons: while math- 
ematically similar they are associated with 
quite different physics. The "particle" dyons 
are time-independent 3-d objects and their 
fermionic zero mode are 3-d normalizable and 
time independent. For the "instanton dyons" 
we need 4-d normalizable zero modes. Index 
theorems associate the latter ones with the 
topological charge Q of the 4-d theory: thus an 
instanton ((5 = 1) consisting of Nc dyons pos- 
sess only one of those, which need to be some- 
how shared between the dyons. Another impor- 
tant technical difference is induced by the fact 
that in QCD-like theories the role of the Higgs 
boson is played by A4 rather than scalars or 
pseudoscalars. As a result, the corresponding 
gamma matrix for Higgs is 70 (rather than 1 or 
75): this makes the interaction with the Higgs 
chirally symmetric. 

Let us generalize the fermionic (-) and 
bosonic (+) boundary conditions to a general 
"anyonic" phase 

V^(/3) = exp{-i(l))ij{0) (22) 

which should be satisfied by a normalizable so- 
lution of the Dirac equation 

p^ = 0, (23) 

containing the gauge field in the hedgehog 
ansatz [51^ 

4" = ^a^JAr^ , (24) 

Al = Ur^ . (25) 

We use the gamma matrix convention 7* = 
cr^ (g) (7% 7^ = cr-^ (g) 1, SO that 75 = 7^727^7^ = 
cr^ 1, and do the calculation for the right 
spinor component. The Dirac equation then 
reads 

(a^)a/3(l^MW(V^i?)f =0. (26) 

where we explicitly wrote the Dirac indices a, /3 
and color indices ^4, 5, and where = (1, ia*). 
Now we ansatz [30l [31] [52 
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We consider the matrix 

VAa 



and ansatz 



(28) 



(29) 



The rule of acting with a color and the spin 
sigma matrices on this object is such that we 
multiply by a color matrix r from the left, and 
if we multiply by a spin matrix a, then we mul- 
tiply from the right, and put a minus sign, i.e. 



aip = r]ea = —rjcre . 
The fermion density is given by 
rt^i = Tr (,?tr,) 



(30) 



(31) 



We now plug the ansatz into ( 26 ) and obtain 
the following two equations 



«2« 



H + 2A 



2 

'H-2A 
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a2 = , 



a2 



(32a) 

ai = . 

(32b) 



where we have assumed ipR (X e-^^^/^, i.e. that 
the Fermion has arbitrary periodicity condition 
in the imaginary time direction. 



(Xiir) — solid, a2(r) - dashed 
1.0 



-0.2 
-0.4 




FIG. 3: The profile of (unnormalized) zero 
mode components ai,2 (the solid and the dashed 
curves) as a function of the distance from the 
dyon. We show four different values of = 
0,0.2v//3,0.4^;//3,0.5v//3,0.55v//3. Note that the 
zero mode delocalizes at = 0.5v/l3 



We solve the Eq. (26) in the Appendix, and 
the result is shown in Fig. [3j However here 



we can easily look at the asymptotic behavior 
of the solution, i.e. when 1-L{r ^ oo) = v and 
A{r ^ oo) = 0, then 



a;(r) + -ai(r) + ^a2(r) = 0, (33) 

4(0 + ^«2(r) + ^ai=0. (34) 

This equation is easily solvable by taking the 
substitution a± = ai ± 0^2 we get 



a± = e 



(35) 



In order for the solution to be normalizable, 
we must have that both a± vanish at infinity. 
This is only possible if < \v\l3/2 



2. The zero mode hopping 

The formulas derived in the previous section 
explain why the zero mode "jumps" from one 
dyon to the other with the change in the pe- 
riodicity condition of the fermions. However, 
there we assumed a static M-type dyon, which 
has the zero mode 



H . 



■101 



(36) 



where all dimensionful quantities are expressed 
in units oi j3 — 1/T. Now it is quite clear that 
(j) G {—v/2^v/2) will preserve the normalizabil- 
ity of the solution. But for the phase in the 
interval v/2 < (/) < tt, one can use the equation 
for the zero mode on the L dyon instead. To 
do so one must first go to the static gauge in 
which V ^ V = 2t: — v^ and then the zero mode 
has good asymptotic behavior 



(M_|^|)^ 



(37) 



Furthermore, because one has to apply a time- 
dependent gauge transformation to reinstate 
^4 — ^ at infinity in the form U = exp(z7rtr -r), 
the fermions (in the fundamental representa- 
tion) gain the desired antiperiodicity. One can 
replace ^ by ^ = tt — ^, or 



(M_|0|), 



(38) 



Now we insert v = 27t — v and (j) = ^ — tt, where 
(j) and V are the (true) holonomy and periodicity 
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of fermions. We assume (j) G (t'/2,7r), so that 
the exponential term becomes 



1^1 



V 

2 



■ TT + ( 



V 

2 



(39) 



and the fermion zero mode becomes normahz- 
able for G {v/2,7r). 

Finahy, we explore the region ^ G (tt, 27r — 
v/2). Then we argue that the zero mode on 
the L dyon is still normalizable. Indeed, the 
exponent now becomes 



Therefore, 



27r - - - 0G (0,7r-'?;/2) . (40) 



-(M/2-|0|)r 



-(|0|-|^|/2)r 
-(27r-|t;|/2-|0|)r 



-|^|/2<0<|7;|/2 



\v\/2<(f)<7V 
7T<(p<27T-\v\/2 



(41) 



C. The fermionic interaction among the 
clusters 

As it has already been mentioned above, on 
general grounds one expects the LL clusters to 
repel each other, as say atoms do, because of 
the Pauli principle. In this section we will show 
how it works using the first "nondiagonal" dia- 
grams in which fermion exchange between such 
clusters takes place. 

The fermion determinant will be of the form 

det^ = |Tn|2|T22r..|r^^|2 

~ ^12^22^21^11^33 • • • ^iViV • • • 5 (42) 

where T^j = — T^^. We can interpret the first 
term as a two-loop diagram, with the fermion 
hopping from one dyon to the antidyon and 
back, for each of the pairs 11,22 etc. The sec- 
ond term is interpreted as a one-loop process, 
in which the fermion is hopping from 1 to 2 
and then from 2 to 2 to I back to 1. The de- 
terminant is the sum of all such terms, with 
the appropriate minus sign to enforce the Fermi 
statistics. Note that the zero resulting from the 
cancellation between hose diagrams means that 
a dyon-antidyon pair will be repelled by another 
dyon- antidyon pair at certain distances. 



It is simple to see that chiral symmetry is 
necessarily restored if the ensemble is made of 
dyonic pairs. Then the determinant is dom- 
inated by the near-diagonal matrix elements 
Tj/, where indices /, 7 go over dyons and anti- 
dyons respectively, which are the closest pairs, 
i.e. Tjj « Tjj, for J ^ I. Then the spectrum 
of the Dirac operator is exactly solvable and is 
given by A/ = zb|Tj j|. Therefore very small 
eigenvalues will be given by very small matrix 
elements Tjj of the dyonic pairs. This matrix 
element is small only if the respective dyons 
are very far away (much further than the range 
of the transition element, i.e. l/v). But since 
the overall configuration is anyway weighted by 
the determinant to the power of the number 
of flavors Nf^ these conflgurations are strongly 
suppressed, and the density of such eigenvalues 
goes to zero at small eigenvalues, implying that, 
by the Banks-Casher relation, chiral symmetry 
is restored. 

An 2A^ X 2A^ matrix of the form 




M = 

has a determinant equal to 
detM = |det 



(43) 



(44) 



which is always positive. 

Let us now consider the fermionic determi- 
nant in the basis of fermionic localized modes 
for 2 dyons and 2 antidyons, labeled with in- 
dices 1, 2 and 1, 2 respectively. 



det^ = |TnT22-Ti2T2i| 



(45) 



As an example, consider a conflguration of 
dyons and antidyons placed on a rectangle of 
dimensions a x b. A little thought will imme- 
diately reveal that if we put them on a square 
such that as we go around we have 1122, the 
determinant vanishes when a = 6, or in other 
words clusters 1 — 1 and 2 — 2 are mutually al- 
ready inflnitely repelling when b = a. However 
we can make them come closer if we orient them 
on the rectangle as 1122, i.e. dyons 1 — I and 
2 — 2 form independent clusters with distance 
^11 = = a and = = y/la^ + b'^). 
Then the repulsion for small b/a will be 



Veff = -ln(det^) - -log[r(a)r'(a) 



(46) 
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where T{rjj) = Tjj. Quite clearly the effective 
potential becomes infinite when 6^0, making 
an effective repulsive core for two dyonic clus- 
ters. 

To discuss this further we introduce the di- 
agrammatic interpretation of the determinant 
viewed as a sum over all fermionic loops. Let 
us view a determinant in some basis of local 
fermionic states. This need not (and in fact 
most certainly is not) be an eigenbasis of the 
Dirac operator. The basis vectors we denote as 
tpn^ which are localized at Xn- Since this basis 
is not an eigenbasis of the Dirac operator 
we have that 



The interpretation of this expansion is then 
straightforward. The determinant can be inter- 
preted as loop diagrams connecting the various 
dyons which carry a zero mode. It is quite ev- 
ident from this diagrammatic expansion that 
every diagram of two loops will have a similar 
diagram with the opposite sign where the two 
loops join via a small channel (see Fig. [4|. 



Jn 



(47) 



where Jn is a spinor resulting from the action 
of the Dirac operator. However we may view 
Jn as a source of our basis states, and interpret 
V^n(^) = / d'^yA{x - y)Jn{y), where A(x - y) 
is the fermionic propagator. Then p taken be- 
tween two states V^rn,n will be 



d'^xd'^yj}^{x)A{x - y) J„(y) . 



(48) 

Therefore we can view the matrix elements 
Tjj as being integrated propagators from one 
source to another. 

The diagrammatic description of the deter- 
minant A in the upper right quadrant is then 
(we assume Nf = 1 in what follows) 



detA 



'■■■T+(all possible perm) 



where the black circle represents the dyon and 
gray the anti-dyon. 

The complex conjugation can be viewed, in- 
stead of fermion going from the dyon to anti- 
dyon, as the opposite propagation of an anti- 
fermion going from an anti-dyon to dyon. The 
pictorial representation of the determinant is 

detp = \ det = det(A) det(A'^) = 





FIG. 4: (Color online)A graphical interpretation 
of the weight in the background of the dyons. Note 
that the relative minus sign will always induce a 
repulsion between a dyon-antidyon pair. 

Let us now think how many pairs can we 
make from N dyons and N anti-dyons. After a 
bit of thought we can see that it is N\. All terms 
like that involve permuting in the above ex- 
pression two positions of two (anti) dyons, and, 
because it requires two exchanges, the sign re- 
mains the same. The determinant will be in- 
tegrated over the moduli space of the dyon- 
antidyon pairs, and so all of these kind of per- 
mutation can be taken to be the same. There- 
fore the fist term will contribute with a factor 
of A^! to the determinant. 

The second term is a bit more tricky. We 
proceed in the following way. Let us consider 
k 4-plets (a 2 dyon-antidyon pairs) which fa- 
cilitate one loop. The number of ways we can 
have one 4-plet is (A'(A^ — 1))^, because we can 
pick two dyons in N{N — 1) ways, and the same 
for anti-dyons. For k such 4-plets we have the 
expression 



(49) 



_ [N{N-l){N-2)...{N-2k + ll 

(50) 

where the k\ factor is present to compensate 
for all possible interchanges of all k 4-plets. 
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The rest of dyons and antidyons can be made 
into pairs, and since there are N — 2k leftover 
(anti)dyons we can do this in [N — 2k)\ ways. 
Combining this with the above factor we get 



2^{N -2k)\k\ 



(51) 



Now we consider the integrals over such ma- 
trix elements. Generally we will have that 
for the arrangement of pairs the integral over 
the moduli space (assuming flat moduli space 
metric - an assumption justified only in dilute 
phase) will be given by 



iV! Co 



V 



N 



(52) 



where cq is a constant which depends on a 
particular form of the matrix element Tjj. 
In other words we have written the integral 
j{T{r))'^d^r = co/m^, making it explicit that 
the effective volume of the integration measure 
is given by 1/m^, i.e. integrating over pairs will 
introduce a volume given by the range of their 
matrix elements Tj^i ~ g-mr^ overall 
volume corresponding to the integration over 
the center of mass of each pair. Notice that we 
also put in the factor A/"!, which is an overall 
degeneracy of the integral. 

In the case of k 4-plets and N — 2k pairs, we 
have 



V 



2^{N -2k)\k\ \^^'m 



V 



N-2k 



(53) 



The constant Ci appears in the quadrupole in- 
tegral over a loop which includes 4- ( ant i) dyons 
(2 dyons and 2 antidyons). The effective vol- 
ume is now 1/m^, with a single, overall, volume 
factor. Therefore the partition function can be 
approximated as 



V 



NN/2 

E 

k=0 



where n = N/V. Rewriting 



(-l)'=Ar! / cin 
2'=(A^ - 2A;)!fc! \4mM 
(54) 



N\{co 



V 



N N/2 



(-l)'=iV! 



^ (A^-2A;)!fc! \N 



A 



(55) 



where the factor 

A = {cl/2c0^)n/m^ 



(56) 



The coefficient cl/cO^ is just a numerical factor, 
and it depends on how fast the matrix element 
falls of with distance, but it does not depend 
on overall coefficient in front of the transition 
element. 

The sum above is can be computed by using 
the identity 



N/2 



m 



k\{N -2k)\ 



{2x) 



N-2k 




Then the partition function becomes 

The approximation can only be valid if A is 
small, therefore 1/A is large. We can employ 
an asymptotic form of Hermite polynomial for 
large N in the following form ^48j 



2iV/2-3/4y^ 

{7TN)i ^/s'mh(|) 



(58) 

,(iV/2+l/4)(20-sinh2(/>) ^^g^^ 



for X = \/2N + 1 cosh (j). Another asymptotic 
series assuming x = \/2N + 1 cos (j) leads to os- 
cillatory asymptotics, which is clearly a good 
indicator that our approximation of just in- 
cluding 4-plet diagrams is invalidated, and that 
higher order diagrams become important. For 
such an estimate the chiral symmetry will be 
restored for A < 1/8, or n/m^ < Cq/2ci, where 
n is the density of one species of dyons. 
^ Finally going back to ( [54| ) for a moment, 
we see that if we look for the quantity (k) = 
Ad{\n Z)/dA^ each coefficient will have a factor 
of k in front. (Note that this is not and aver- 
age number of 4-plets: each configuration has 
arbitrarily many 4-plets, 6-plets, etc.)In Fig. Is] 
we show (k) /N as a function of parameter A. 
Notice the abrupt change as we approach the 
critical value A = 1/8, which we take as an 
indication of chiral symmetry breaking. 
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FI G. 5 : (k) /N as a, function of parameter A defined 
in (l56l . 



D. 



Bosonic One-Loop Interactions and 
Electric Screening 



Tfie basic pfiysics of tfie electric screening can 
be explained most simply following the original 
derivation by one of us [33] (in the Coulomb 
gauge). If some object possessing a nonzero A4 
is immersed into a quark-gluon plasma, those 
quanta from the heat bath are scattered on 
it. The simplest diagram comes from the quar- 
tic term in the gauge Lagrangian, ~ g'^A^A\ 
which couples the heat bath gluons directly to 
square of ^4, but there are also other diagrams 
contributing to the forward scattering ampli- 
tude. The result was the expression for the 
QCD Debye mass [53^ 



D 



g'T\N,/3^Nf/6) 



(60) 



for massless quarks and gluons: incorporation 
of their effective masses is straightforward. 

The next relevant paper is by Pisarski and 
Yaffe (PY) [3T who calculated the one-loop ac- 
tion of the calorons (the finite-T instantons). 
Its main part is the following correction to the 
instanton action 



2^2 



SSpY 



D 



(61) 



where p is the instanton radius. The first factor 
in this expression comes from the (4d) dipole 
moment of the instanton, and the second from 
the forward scattering amplitude of the ther- 
mal plasma quanta on it, for derivation see [35]. 
This term is only present in the plasma phase. 



at T > Tc, as only in this case there exist ther- 
mal quarks and gluons undergoing this scatter- 
ing. 

Going forward to calorons at nonzero holon- 
omy, a corresponding one-loop effective ac- 
tion has been computed by Diakonov, Gro- 
mov, Petrov and Slizovskiy (DGPS) [36]. The 
caloron is now a superposition of the M and 
L dyons, separated by distance tml, and the 
basic expression from which the effect comes is 
the following integral 



J 



1 



-f = Atttml + . 



(62) 

where vl^vm are distances from the dyon cen- 
ters to the observation point r. 

(The dots stand for corrections due to a fi- 
nite dyon size: the Coulombic A4 is true only 
at large distance. Note also that at large r 
the integral converges because the integrand is 
~ t'^m/'^^' This term comes again from the 
quartic term in the action, in which two gauge 
potentials are the A4 of the instanton and two 
others belong to the thermal gluons.) 

Thus, the electric screening effectively gener- 
ates the confinement of two dyons, with a po- 
tential linearly depending on the LM separa- 
tion: 



V 



LM 



27rMl 



(63) 



At the zero holonomy this result matches the 
PY answer because of the "instanton size rela- 
tion" 



np T = tml 



(64) 



which, so to say, relates the 4-d dipole of the in- 
stanton to the 3-d dipole of the dyon pair. The 
second factor is still the same thermal integral 
producing T^. One obvious effect of fermions 
is that a generalization of the DGPS result to 
theory with the generalized Debye factor as in 
(60). As a result, the electric screening effect 



ensures LM "binding" into a finite-size object 
with the (inverse) size 



T{Nj3^Nf/6) (65) 



Note also, that for LL or MM pairs with 
the same electric charge, there will be a plus 
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in the integral (62) above and thus the effect 



becomes repulsive and the integral diverges: it 
needs to be regulated by some opposite charges. 
For molecules consisting of all 4 (L, M, Z, M) 
dyons, to be discussed shortly, the screening 
potential is 



Ml 



<A\> 



LMLM 



(66) 



in which the ^4 now contains all 4 Coulomb 
contributions. This integral is of course con- 
vergent because of total the zero charge of the 
molecule. 

Let us give an example of the electric multi- 
dyon screening potential created in this case. 
We will later see that the direct fermionic in- 
teraction binds LL pairs stronger than the LM 
interaction. Therefore, for simplicity one can 
ignore the LL cluster size and put them at the 
same point, the origin. Another simplification 
appears if one puts M, M and LL on one line. 



The integral ( 62 ) changes to 



1 



1 



(67) 



The corresponding potential is shown in FigjSj 
As one can see, like for DGPS case, the poten- 
tial consists of linear segments, but is now de- 
formed away from the companion dyon. (Note, 
that it is not due to their Coulomb repulsion, 
which is also there but will be discussed in the 
next section.) 



III. STATISTICAL MECHANICS OF 
DYONS 

A. Statistical Mechanics of a Single 
Dyonic Molecule 

The partition function for an instanton- 
antiinstanton molecule can be recovered us- 
ing known elements for each of the ingredi- 
ents. Let us start with dimensional consider- 
ations, valid at high enough temperatures. If 
the fermions are all massless, the overall power 
of the T-dependence of total molecular density 
can be determined from the known power of the 
^QCD in the instanton-antiinstanton measures. 




FIG. 6: (Color online) The integral (67) propor- 
tional to the potential for the M (or M) dyon cre- 
ated by the electric screening as a function of its 
position. The charge-2 LL pair is assumed to be 
in the origin, and the companion dyon is placed at 
the point 1. 



namely 



'^mol 



logZ^ 



mol 
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(68) 



Typically the power in the second factor is large 
and negative, so this density rapidly decreases 
with T. (It is so except near the boundary 
of the asymptotic freedom domain where that 
power is getting small: we will not discuss this 
region.) 

Translation from the dyonic to instanton lan- 
guage at the level of the moduli metric and par- 
tition function has been studied for self-dual 
and anti-self- dual sectors by Diakonov et al 
[36]. Their expression, in the SU(2) case for 
pure gauge theory is 



dZLM = d^rLd^rMT''2nC 



(69) 



47rT 



where tlm = \^m — r/vfl, numerical constant 
C = 1.0314..., and the scale parameter Apy 
is for the Pauli-Villars regularization [54] The 
factor 



FD(r)= (1 



vvr \ 
2^) ^ 



+ vr) ^""^ (1 + vr)^^ 



is the correction appearing due to non-zero 
holonomy. If the holonomy ^' = or antiholon- 
omy ^ = 0, in the expression above Fjj = 1 
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modes as well as the nonzero mode part 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.! 

FIG. 7: Histograms of the distributions of the 
distance between the LM dyons (stars) and LL 
(boxes) dyons, for A^c = 2, A^/ = 4 molecule, in the 
units of the Matsubara time 1/T. The fermionic 
mass is taken to be m/ = T and the holonomy 
ly = 0.1. 



and it reduces to the well-known caloron mea- 
sure, using the relation (64). In the limit of 



large dyon separation one may keep only the 
r-terms: note that it then becomes fiat and r- 
independent as one would expect. The screen- 
ing potential for LM pair is 



Vscr{r) 



27rr 



.T(l--L)-. 



(70) 



-.T(l--^)+. 



We have excluded one more factor in the parti- 
tion function of [36l 



exp{—V 



(71) 



which does not depend on the calorons/dyons 
and is just a one- loop contribution to the prob- 
ability to have the holonomy v in the ensemble: 
it certainly should not be repeated twice. 

The same expression can be repeated for 
the for the Z, M pair, and then the com- 
bined into the 4-particle partition function for 
a "molecule". Since, unlike [36 , we are inter- 
ested in the theories with fermions, we intro- 
duced extra factors which include that for zero 



dZmoi = dZ^MdZ 



m2 + \TiA{rLL) 



A2 



Nf 



(72) 



C{Nf) 



Nf/6 



2^2 



As discussed in the preceding section, Vgcreening 
is defined by the 4-particle expression for A4 
integrated over the volume. If one of the par- 
ticle is going to large distances, the expression 
reduces to a dipole and return the linear con- 
finement result, preventing "ionization" of a 
molecule. 

The bracket in the power of Nf/6 comes from 
the nonzero mode part of the fermionic deter- 
minant calculated by 't Hooft. The power of 
A in it corresponds to correct beta function of 
the theory with Nf fermions, and is therefore 
fixed. Its dimension should be compensated by 
some parameters with the dimension of the dis- 
tance, which in the case of a single instanton 
can only be its size p. At finite temperature 
instantons lose 4-d spherical symmetry and an- 
other dimensional parameter - T appears, as 
well as a nonzero holonomy v. Lacking explicit 
evaluation of the nonzero mode determinant, 
we just used the instanton expression and the 
relation irp'^T = r translating its size into the 
current language. At least for small-size instan- 
tons (small r) this should work. As a result we 
get factors r^^f^'^ in the measure, or a repul- 
sive potential ~ Nflog{l/r) trying to dissolve 
the molecule. Recall however that it only sup- 
pose to be true at small r, while at large r the 
one-loop electric screening effects generate an 
attractive potential linear in r which would pre- 
vent it from happening. 

We have introduced here the fermion mass m 
in the Dirac operator, assuming it is the same 
for all flavors, for the normalization reasons [55 . 
The term proportional to the masses is nothing 
else but a square of the independent instan- 
ton and antiinstantons: and since their nor- 
malization have been already determined by 't 
Hooft, the flavor-dependent normalization con- 
stant C{Nf) can be determined for Apy. 

If the fermion masses are set to zero, the 
fermions couple instanton to antiinstanton via 
the overlap matrix element and four integrals 
over dyon positions produce three convergent 
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integrals, while one remaining integral over the 
global position produces one factor of the 
box volume. 

Even for the simplest case of SU{2) color, 
when molecules contain 4 dyons, their position 
space is already 12-dimensional. Therefore we 
used standard Metropolis algorithm to gener- 
ate their statistical distributions. FigjT] shows 
one typical example of the output: in it we 
compare the distances between the LM dyons 
(stars) with that of the LL (boxes). The latter 
are seen to be much tighter placed, forming a 
"nucleus" of the molecule. 



B. Modelling the dyonic ensemble 

1. Three Molecular Models 

As a first step toward the understanding of 
the dyonic ensembles, and their role in chiral 
symmetry breaking/restoration, we had formu- 
lated some simplified models. 

For calculation purposes it is convenient for 
these models to treat the dyon density 



given by 



rid riL 



riM 



(73) 



(which is also the same as the instanton density 
riinst + riantiinstanton) as the basic dimensional 

— 1/3 

quantity, providing the units of length ' . 
Using such length units we put n^^ = 1 for a 
while, and will be expressing other dimensional 
quantities in these units. We will be working 
with traditional periodic boxes of some size L x 
L X with L "large", and thus put into such 
boxes Nd = dyons of each kind. 

For each configuration of these models we 
then calculate the fermionic matrix T^j, and 
calculate its eigenvalues. In this way we get 
part of the Dirac spectrum built on the sub- 
space of the dyon zero modes. Since antiperi- 
odic fermionic zero modes resign on L, Z dyons 
only, fermionic part of the measure ignores the 
M, M dyons. The matrix is thus of the size 
2Nd X 2Nd. For reasons of opposite chirality, 
two quarters of the matrix, when both z,j = L 
or L are zero, so fermionic hopping occurs only 
from a dyon and anti-dyon. 

We assume the matrix element T^j, where is 



-Mvij 



(74) 



where one of the indices counts L dyons, and 
the other counts L dyons. The constant c will 
be left undetermined, whereas the "mass" M is 
given by 



M = v/2 = (27r - v)/2 . 



(75) 



The form of the matrix element is not de- 
rived, but it is postulated as expected form for 
the matrix element of the Dirac operator in the 
zero mode basis. We introduce a regulator at 
r = 0, so as to make deal with smooth dis- 
tributions. However there will be an natural 
cutoff for how close the dyon and antidyon can 
get before they are part of a perturbative vac- 
uum, which is roughly the size of their cores 
1/M ~ 1/v. Strictly speaking the definition of 
the distance at which the dyon-antidyon pair is 
irrelevant is defined as the distance at which it 
no longer supports a localized fermionic mode. 

We proceed by three models: 

1. The Random Gas Model 

2. The Random Molecular Model 

3. The Reweighted Molecular Model 

The simplest model-I is that of the "Random 
Gas Model", in which all correlations between 
the dyons are ignored and they are placed ran- 
domly. The only parameter of the model is the 
fermion mass M entering the matrix T^j, to be 

1 /3 

expressed in units of nj . (In reality, both the 
dyon density and the holonomy, defining M, 
will be function of the temperature T, but we 
prefer to study our models in their parameter 
space before mapping some of the results to lat- 
tice data, see below.) 

In Figs. [8] and [9] we show the results of such 
calculation. We use the box of the size 6^ and 
thus 216 dyons of each kind, and a range of 
fermion masses as indicated in the figure cap- 
tion. The characteristic feature of the "Ran- 
dom Gas Model" is a large peak near eigen- 
values A ^ 0. Since the density of quasizero 
eigenvalues is proportional to the quark con- 
densate (Casher-Banks theorem), we conclude 
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FIG. 8: (Color online) The Dirac eigenvalue spec- 
trum for the random dyon gas 



that this model provides large or "enhanced" 
chiral symmetry breaking. 

Our 2^^^ model is the "Random Molecules 
Model", in which we include pair correla- 
tions between L, Z dyons. As we discussed 
above, we expect significant attraction between 
those, of classical (Higgs-related) origin, as well 
as the fermion-induced dion-antydyon confine- 
ment. We focus here on the fermionic interac- 
tion. As the number of fermionic zero modes 
grows, proportionally to the number of flavors 
Nf^ we expect that at large enough Nf the 
molecule mean size Rm decreases as ~ 1/^/- 
We model the vacuum as being composed out 
of random molecules. The following distribu- 
tion of the size of molecules is used: 



Dmoi{r) = Nr^ 



-Mr 



vTTMr 



2Nf 



(76) 



where is due to the measure of the dyon- 
antidyon coordinates, and N is the normaliza- 
tion constant. The above form is inspired by 
a weight (detT^j)^^, for a dilute molecule en- 
semble. Average molecular size will roughly be 
given by Rm = l/{Nfm). 

At this stage we have ignored any interaction 
between the molecules, placing them randomly 
with random orientations. 

The model has two parameters, the holon- 
omy "mass" M and the number of flavors Nf. 



In Fig. 10 we show the Dirac spectrum for sev- 
eral values of Nf^ as a function of M and in 
Fig. [iT] the lowest eigenvalue and chiral con- 



The 3^^ model is a modification of the second 
by reweighting the configurations with the de- 
terminant {det Ti j)'^^ ^ / (Yl-Dmoii'^i^))^ where 
ru is the distance between closest neighbors. 
The result is shown in Fig. [inland in Fig. [13] 



2. Mapping the models to finite T QCD 

Our three models were defined in such a way, 
that each step has introduced one new param- 
eter: with the dyon density it brings the total 
number of parameters to four. Yet a QCD-like 
theories with massless fermions have only two 
parameters, the temperature T and A. Thus 
only a 2-parameter subspace of our (up-to) four 
parameter model space can be compared to re- 
ality. 

Now is the time to map those parameters. 
Like lattice practitioners do, it is thus natural 
to measure all dimensional quantities in units 
of T. The dimensionless dyon density Ud/T^ is 
one of the key parameters. It has not yet been 
measured on the lattice, but it can be. While 
it is the same as the density of instantons, it is 
not given by the topological susceptibility x{T), 
as neutral molecules contribute to the former 
but not the latter. Semiclassical theory tells 
us the large-T asymptotical dependence on T, 
expression (68). For qualitative estimates one 



may use it normalized to its value at Tc . The 
factor in front of the power oiTc/T depends on 
the coupling, in a particular definition used by 
'tHooft, and the fermionic factors its value is 
0(1) for physical QCD or Nf = 2, which one 
can use for absolute normalization. 

The fermion mass M/T ^ ^/Nf and thus, 
keeping the coupling fixed while increasing Nf^ 
one finds that the cluster size is reducing and 
thus we are going into a regime of more dilute 
gas. However, if one wants to follow the line 
of "fixed eigenvalue spectrum" and/or fixed < 
>, one needs to keep the same diluteness of 
the molecular model, or keep constant 



RmTid = const , 



(77) 



densate results. Note that here there is explicit 
dependence of chiral condensate on holonomy. 



where R^ ~ l/{NfM). Thus the dyon density 
should grow as Ar|, e.g. from Nf = 2 to 12 
increase by a factor 216. 

The only way it can be achieved is by a shift 
into the stronger coupling! A crude estimate 
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FIG. 9: (Color online) The eigenvalue gap and the quark condensate for the random gas model, as a 
function of the "holonomy mass" M. 
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FIG. 10: (Color online) The Dirac eigenvalue spectrum. The horizontal axis is the eigenvalue of the Dirac 
operator A expressed in units of c defined in the eq. ( |74[ ), and vertical axis is the probability density dP/X. 
The plots are for Nf = 2 and M = (tt/G . . . 9) x 7r/6rr^, where n is the density of L dyons. Note that the 
chiral symmetry is restored as a function of M which is connected to holonomy as M = ^;/2 = (27r — v)/2. 



ignoring preexponent gives a shift of 

(78) 

This qualitatively explains why chiral restora- 
tion line < i/ji/j >= derived in lattice stud- 
ies (see Figj2| dramatically shifts into stronger 
coupling. 



Unfortunately our attempts to do it quanti- 
tatively failed, for the following reason. The 
coupling g in the semiclassical expressions and 
on the lattice (such as Figj2| are defined in dif- 
ferent schemes, with the scales and A^^f 
Perturbatively they only differ by a calculable 
factor, but as their ratio happen to be large, 
this relation is not very useful in practice. For 
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FIG. 11: (Color online) The smallest eigenvalue and the chiral condensate for the Random Molecule model, 
as a function of the holonomy mass combined with the dyon density, for Nf — 2. 
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FIG. 12: (Color online) The Dirac eigenvalue spectrum for a Random Molecule Ensemble - light (blue) 
and Reweighted Random Molecule Ensemble - dark (red). The horizontal axis is the eigenvalue of the 
Dirac operator A expressed in units of c defined in the eq. (74), and vertical axis is the probability density 



dP/\. The plots are for ATy = 2 and M = (30,25,20, 15, 10,5) x 7r/6n^/^, where n is the density of L 
dyons. The reweighting becomes unreliable in the last plot, and only one configuration dominates. 



example, for QCD or Nf = 3 theory, the former 
is about 300 MeV and the latter about 5 MeV. 
The instanton density includes huge factor 



exp 



8^ 

n2 



(79) 

Obviously, in view of such a huge factor, any 
small deviation from the two-loop beta function 
used would result in huge uncertainties which 
make any numerical comparison of the semi- 
classical expressions and lattice bare coupling 
meaningless. 

The tests of this explanation however can 
still be made using lattice data. The most 
straightforward one would be to measure the 



dyon molecule density and size and see if the 
relation (77) holds. To do so one can e.g. use 



Dirac eigenstates in certain interval of A which 
can identify a dyon-antidyon cluster. 



3. Structure of strongly coupled ^'cluster liquid" 

With large number of fermions LL are 
strongly coupled into a charge -2 well local- 
ized objects, compensated by the negatively 
charged M and M dyons which are more ho- 
mogeneously distributed. 

The LL "nuclei" , with the fermions attracted 
to them, form mutually repulsive "atoms" . The 
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FIG. 13: (Color online) The smallest eigenvalue and the chiral condensate for the Reweighted Random 
Molecule model, as a function of holonomy for Nf = 2 and M = 7r/6...30 x 7r/6. Note that reweighting 
becomes unreliable at around M = 15 x 7r/6, and that we show this result only to demonstrate the trend 
that chiral symmetry persists to lower values of the holonomy mass M than for the nonreweighted case. 



question is what arrangements those should 
have to get the lowest energy. One obvious idea 
is those with the best packing in 3d, namely 
the face-center-cubic (fee) or the hexagonal- 
close-packing (hep). Selecting between those 
one may follow the guidance given by ordinary 
atoms which are neutral and spherically sym- 
metric by themselves. High-density solid He"^ is 
of the hep structure, and perhaps that would be 
the approximate local symmetry of our strongly 
correlated LL liquid. If so, each of them has 12 
nearest neighbors, organized in two hexagons 
at two planes above and below the cluster. 

While large number of correlated neighbors 
reinforces the correlations, simple estimates 
show that the parameter F (the average inter- 
action potential divided by T) does not reach 
the critical value needed for the ensemble to get 
solidified. Thus, with the parameter range at 
hand we expect a strongly coupled liquid, with 
smaller number of well-correlated neighbors but 
with their locations still correspond to those in 
the crystal. Explicit statistical simulations of it 
are possible, but are deferred for further stud- 
ies. 



C. Statistical mechanics driven by bosonic 
moduli space metric effects 

Let us now discuss the opposite limit of zero 
Nf = 0, in which there are no LL clustering. 



Let us further imagine that for some reason 
one can decouple the two sectors, dyons and 
antidyons, and discuss what kind of a system 
would be created by "Diakonov's determinantal 
forces". Assuming that the antiselfdual sector 
do not exist, let us thus focus on the M, L sec- 
tor. Since the electric and magnetic charges in 
it have the same sign, one may in this section 
simply call them + and — dyons. 

At large distances the forces between them 
are Coulombic, and one may think that the 
local crystal correlations those generate is a 
simple cubic crystal of alternating charges, like 
e.g. the usual salt NaCl. Any charge is thus 
strongly correlated with 6 nearest neighbors. 

Classical Coulombic systems, are well-known 
to be unstable against charges falling on each 
other. (Of course for real ions electron repul- 
sion solves this issue, stabilizing the salts.) We 
thus studied the following question: can the Di- 
akonov's determinantal forces stabilize a cubic 
crystal? 

We use the moduli space metric Gij for the 
selfdual sector as suggested in [2] to calculate 
the effective potential of a crystal configuration 
of L and M dyons, with a lattice spacing a, for 
a displacement of a single dyon somewhere in a 
center. Effective potential is the log of the mea- 
sure, Vd = — In(detG). For a purely Coulom- 
bic crystal, the crystal potential as a function 
of a displacement Ax has infinite Coulombic 
dips as displacement Ax approaches ±a, which 
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means that for purely Coulombic interactions 
the alternating charges will fall on each other. 
However the effective potential Vd contains re- 
pulsion and, as shown in [l4j this leads to a 
pronounced minimum at x = for sufficiently 
small lattice spacings a (high density). On 
closer inspection, however, there is always a 
small, but clear minimum. As expected, the 
Coulombic dips at ±a still persist. This diver- 
gence -corresponding to small-size instantons 
and the factor in the measure - is known 
to be removed by quantum fluctuations, which 
produce a stronger factor pii^c/3-2Arj/3 ^j^^g 
outside of classical approximation the problem 
is well defined. 




FIG. 14: (color online) The effective potential 
Vd = — In(detG'), where G is the Diakonov's de- 
terminant, as a function of the displacement Ax 
of a single dyon in the center of a cubic 6x6x6 
crystal, in the direction of the adjacent dyon. The 
graphs have been rescaled for better comparison as 
follows: (blue) solid 684^^, a = 0.1; (red) dash 818 
Vd,a = 0.25; (brown) short dash 917Vd,a = .5; 
(green) dash-dotted 967Vd,a = 0.75; (blue) long- 
dash 1000 Vd, a = 1. The units of Ax, a is the Mat- 
subara time. Note that as dyonic density increases 
by a factor 10*^, this one-loop bosonic interaction 
creates a significant minimum at AX = 0, stabiliz- 
ing the cubic structure. 

On the other hand, as discussed in [2^ it 
seems that from the point of view of the far field 
metric, the antiselfdual sector behaves similarly 
to the selfdual one, and the interaction between 
L and L is similar to the interaction between L 
and L (i.e. repulsive), while that of L and M 
is attractive. Therefore we have three possible 



structures depicted in Fig 15 



fermions correlate the L and Z, then they will 
repel the other pairs, therefore making the 
hexagonal crystal, as we discussed above. For 
zero Nf they can be either strongly correlated 
(b) or form an alternating crystal (c). Lattice 
practitioners often introduce the so called "va- 
lence" quarks, which are not in the measure but 
are only used for a diagnostic purposes of the 
"quenched" {Nf = 0) theory. Dirac eigenvalue 
spectrum and chiral properties revealed by such 
studies can be computed and compared to the 
lattice data. 

In a standard way we model the Dirac ma- 
trix by a "hopping" matrix ( |2|3[ ), with the ma- 
trix elements being some function of a distance 
between dyon-antidyon /(r). We expect these 
functions to be exponential in distance r = \ f\ 
at large r and some constant we call /(O) at 
small r. 

In the first approximation we can consider 
only next-neighbor matrix elements, and use 
the cubic structure (b). Then we can write the 
upper-right (or the lower-left) part of the Dirac 
operator matrix, in the triple-index notations 
with n, m, / being positions in units of a along 
the three spatial coordinates in the cubic lattice 
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Upon standard diagonalization by transfor- 
mation to the dual lattice momentum states 



ifc) 



-{n,m,n) 



(81) 



we obtain that the spectrum is given by 



iy{k) = /(0)+/(a) (cos h + cos(A:2) + cos{ks)) , 

(82) 

where ki^2,3 go from (0, 27r), and the elementary 
number of states is given by standard dN = 
Vd^k/{27r)^. The density of states is 



dN 
'dX 



V J d^k S{X — cos ki — cos /c2 — cos ks 



(83) 



For nonzero number of quark flavors Nf 



where we have used the shifted eigenvalues 

A = (^^ - /(0))//(a) (84) 
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The spectrum can be integrated to yield 

1 



tion for chiral symmetry breaking is therefore 



dN 



I 



dkidk2 



smks{ki,k2)\ 



(85) 



where ks = arccos(A — cos/ci — cos/c2), 
and the region of integration is such that 
— cos /ci — cos A:2| < 1. Numerical integra- 
tion yield the curve shown in Fig. [16] 
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FIG. 16: (Color online) The plot of the density of 
states of the cubic lattice with the next-neighbor 
interactions. 



We see that the density of states form a 
band with a sharp boundary, it goes to zero at 
|A| > 3. For scenario (b) this shape will appear 
in the spectrum centered around ±|/(0)|, and, 
for vanishing /(a), will be delta function-like. 
Each of those morphs into a shape of Fig jlH] in 
a type-(b) crystal as /(a) increases. If they are 
separated by more than the width of the peak, 
the chiral symmetry is not broken. The condi- 



/(«) 



/(O) 



> 



1 



(86) 



Alternative structure (c) can be motivated as 
follows. As mentioned earlier, the long range 
interactions between Ls and Ms, become the 
same regardless of them being dyons or an- 
tidyons. This means that the cubic crystal will 
have Ls and Zs located at positions of " + " ion 
and M and M at positions of "-" ion. In this 
case, the spectrum of the Dirac operator is 

v = 2f{aV2){cos ki + cos /c2 + cos k^, + cos k/^) , 

(87) 

where we have approximated that, on average, 
each L has 4 nearest Zs. Similarly as before we 
get that 

^(. = 0) = 2/(av/2) /^^l^ 
dv^ ^ ^ J (27r)4 |sinfc4 



where k^ = arccos(cos ki + 
cos ki + cos /c2 + cos ks < 1 



cos k2 - 



(88) 
cos/ca), and 



D. Note on confinement of the cubic 
crystal 

Here we show that the Polyakov loop of a 
crystal configuration is indeed zero. We have 
that 

trL(x) =tre<5+^^^^)"^ = 2 cos Q + F(a^)) , 

(89) 

where we used the gauge-combed gauge, and 
where V{x) is some potential which goes to 
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— 'u/2 at the position of M and M dyons and to 
— 'u/2 + 7r at the position of L and L dyons. Us- 
ing the identity that cos (a + /3) = cos a cos {3 — 
sin a sin /3, we have that 

V V 

tr L{x) = 2 cos - cos(V(ic))+2 sin - sin{V{x)) . 

(90) 

The above expression has to be integrated over 
ah possible crystal orientations and positions. 
This is equivalent, though, to integrating over 
X in the region of one crystal unit, i.e. 

{ti L{x)) = 2^ 1 d^x [2cos^cos{V{x)) 

CL J ^3 Z 

+2sin^sin(V(ic))] (91) 

However, setting v = (maximally nontriv- 
ial holonomy), we see that the first term van- 
ishes because cos(7r/2) = 0, and the second 
term vanishes because V{x) is alternating be- 
tween L and M dyons in sign, and, since sin is 
an odd function, this term too vanishes. There- 
fore the Polyakov loop averaged over crystal 
configurations vanishes, and it depends explic- 
itly on the holonomy being nontrivial. 

Finally let us note that such "confinement- 
in-average" phenomenon has its predecessors, 
resembling a finite-density holographic model 
of densely packed baryons ^ . In it baryons 
are modeled by instantons, which also undergo 
transition into a "dyon phase" in which they 
restore chiral symmetry in average. 

IV. SUMMARY AND DISCUSSION 

In this work we have done qualitative study 
of the interactions of (anti)selfdual dyons, as 
well as some study of the statistical ensembles 
following from them. We emphasized in partic- 
ularly the importance of the dyon-antidyon in- 
teraction, both classical (bosonic) and fermion- 
induced ones, instead of the interactions inside 
the selfdual and antiselfdual separate sectors. 
We specifically were focusing on the SU(2) the- 
ory and the temperatures above and near the 
chiral restoration phase transition. 

We already summarized the overall picture 
resulting from this study at the very beginning 
of the paper, so let us just outline the elements 
of the picture which we believe explain certain 



lattice observations. We will also emphasize 
what lattice practitioners can do to test our 
predictions further. 

1. Existence of the "topological molecules" 

At high T the instanton/dyon density is 
small due to electric screening, the topological 
object must make "molecules" , neutral in terms 
of all three charges involved: topological, elec- 
tric and magnetic. In the SU(2) theory they 
contain all 4 types of the dyons. We predict 
a peculiar distribution, with "nucleus" of LL 
cluster at the center, with M and M at the pe- 
riphery, see Figjl] This effect is predicted to get 
more pronounced with the increasing number of 
quark fiavors Nf. It is important to check it on 
the lattice, perhaps by identifying the fermionic 
sates "at the gap" in ensembles with varying 
Nf. Its generalization to any number of colors 
is straightforward. 

2. The critical line versus the fermion numbers 

Model simulations with such clusters, ran- 
dom or interacting with each other, predict cer- 
tain distinct shapes of the Dirac eigenspectrum 
at small A. We have in particularly found at 
which diluteness of the "cluster gas" one gets 
particular value of the chiral condensate, the 
chiral symmetry restoration or certain size of 
the "gaps" . These eigenvalue spectra and "lines 
of constant condensate" can be compared with 
the lattice ones, in order to see if chiral break- 
ing does happen in the "molecular gas" regime, 
or at a denser regime. 

The increasing number of fermions leads to 
stronger fermion-induced interactions, binding 
the LL pairs into tighter clusters. In terms of 
our molecular gas model, it gets much more 
dilute, unless the overall density of the dyons 
is significantly increased. This can only be 
achieved by going into stronger coupling do- 
main, which reduces the dyon masses and inter- 
actions. That is why the critical lines in Figj2] 
go downward with increasing fiavors. 

We have also explained why there is a quali- 
tative difference between the fundamental and 
adjoint fermions. While the former have zero 



27 



modes only with one (heavier) dyons LL, the 
adjoint have zero modes for ah dyons, includ- 
ing lighter M, M (for SU(2)). The latter are 
much less correlated, thus their chiral restora- 
tion temperature is much higher. 



3. Chiral splittings of hadronic masses versus 
Nf 

As this parameter is now becoming measur- 
able on the lattice, with the progress in com- 
puter/lattice technology, it is perhaps time to 
map it more consistently, and also think again 
about the physics it reveals. 

Rapid decrease of the chiral condensate scale 
around Nf ~ 4 has been in fact predicted 
by the instanton liquid simulations |7j long 
ago. The reason for that has been a "dip" 
in the eigenvalue spectrum developed due to 
the "molecule" formation [56 . This very phe- 
nomenon is in fact central in our current study. 
What has not been predicted in 1990's was a 
significant shift to the stronger coupling and 
drastic increase of the overall instanton/dyon 
density (or, in alternative language, a rapid de- 
crease of the hadronic scale), which makes a 
very small quark condensate relevant. It would 
be important to study transition to "molecular" 
topological structure in lattice simulation. 

The dependence of the quark condensate on 
the density of "molecules" we found in our cal- 
culations is interesting. As seen in FigjlOj at 
low density there is a minimum between two 
"molecular peaks" , but at some diluteness there 
is a sudden appearance of a small peak inside 
this minimum. This implies a sudden jump 
in the quark condensate value, in a small in- 
terval of parameters. It is more pronounced 
than many crossover phase transitions: and 
thus we may call this phenomenon a phase tran- 
sition in the dyonic ensemble, from "atomic" 
to a "molecular" state, at Nf > 3. Using 
quark masses as interpolating parameters be- 
tween Nf = 3, 4, 5, lattice practitioners can see 
if this change of Dirac spectrum is also occur- 
ring in lattice simulations as well. 



4' Dependence of the Dirac eigenvalues on the 
holonomy value and phase 

The non-zero holonomy provides Higgsing, 
the breaking of the color group, and thus it nat- 
urally explains different fermionic masses which 
appear in the Tij hopping amplitudes. 

Since the rest of the paper has been for SU(2) 
gauge group, let us discuss in a bit more de- 
tail the two SU(3) options shown in Fig 1 of 
[22 , namely the real < P > sector as well 
as the one with the phase of P being 27r/3. 
Generic holonomy in SU(3) is described by 3 
parameters /ii,/i2,M3 subject to one condition 
(/ii + /i2 + = O(mo(i(l)). If one imposes 
one more condition, such as the fixed phase of 
< P >, there is only one free parameter left. 
If the phase is zero, < P > is real, the one- 
parameter family of possible holonomies is 

Ml = 0; /i2 = 1/2 - /is = 1/2 - 5 (92) 

If the phase is 27r/3 this solution is simply ro- 
tated additively to 

Ml = 1/3; /i2 = l/3 + l/2-(5;/i3 = 1/3 + 1/2-^ 

(93) 

The masses of monopoles are determined by the 
differences, which are the same in both cases 

i/i = z/3 = 1/2 + ^,^^2 =2(5 (94) 

since for gluonic observable the two sectors are 
identical by Z^^ symmetry. 

However the fundamental fermions notice the 
difference, as their masses are given by /i^, not 
Vi. Furthermore, which dyon gets the zero 
mode depends on the phase parameter z in 
fermionic periodicity condition: the rule is that 
it is the one in which v sector z resides on the 
circle. The antiperiodic fermions {z = 1/2) 
picks up the second type of dyons V2 = 2(5 in the 
real sector and the first one z^i = 1/2 + (5 in the 
one with the phase 27r/3. The lowest mass of 
the fermion is 27rT(5 in the former case, while in 
the complex ones it is (27rT)mm(l/6, 1/3 — (5). 

Of course, in lattice subsector with the fixed 
phase, the modulus still has some average and 
the distribution, determined by the effective po- 
tential of < P > at a give T, which is known 
if the lattice simulation is made. We however 
dont know the values: let us say take some 
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generic value between and 1 as a guess: a 
half 

I < P > I = 1/2 = (l-2co5(27r(5))/3 0.29 

(95) 

If so, the fermion masses for the two sectors are 
m/(27rT) = 0.29 and 0.04, respectively. Such 
large mass difference explains why the partici- 
pation ratios (roughly, the fraction of the box 
volume occupied by a mode) are so different: 
while in the real sector the Dirac modes occupy 
only about 1 percent of the box, in the complex 
sectors nearly the whole box is occupied.) 



5. Dirac eigenstates "at the gap". 

The objects found via this method by Bruck- 
mann et al ^S] consistent with being made 
of L and L dyons. Apparently thy must be neu- 
tral under the topological charge because they 
do not contribute to topological susceptibility. 
The clusters we propose in this work include 
both L and L dyons, with zero total topology, 
yet still able to support the fermionic localized 
states. Further lattice studies of the abelian 
projected electric and magnetic fields correlated 
with those objects would further clarify their 
origins. 



6. Dependence on fermionic periodicity 
conditions 

In several works [20l [21] effects of the tempo- 
ral boundary conditions on the chiral conden- 
sate on quenched ensembles was explored. In 
addition to the restoration of the chiral sym- 
metry for the physical, antiperiodic, boundary 
conditions, an increase in the chiral conden- 
sate for the periodic condition above Tc was 
observed. 

In the case of the adjoint fermions the drop 
of chiral condensate happens for both the peri- 
odic and antiperiodic sectors. There is however 
a qualitative difference: the drop is slower for 
the periodic case, differs in shape, and was not 
traced to reach zero. A qualitative difference 
is certainly expected on the basis of Ref. [40 
: the zero modes for periodic and antiperiodic 
boundary conditions are drastically different. 



For periodic case, all zero modes are democrat- 
ically distributed, whereas for the antiperiodic 
boundary conditions they all fall onto the heav- 
ies one (the L dyon). The case of antiperiodic 
boundary condition then restore chiral symme- 
try in the way similar to fundamental quarks, 
by condensation into LL clusters. The case 
of periodic adjoint fermions is quite different, 
it has "democratic" distribution of zero modes 
over all dyons of the instanton. This makes 
restoration of the chiral symmetry much more 
difficult, as now correlating the "heaviest" LL 
is insufficient and also "light" MM pairs should 
form. A more quantitative analysis of the ad- 
joint fermion case is postponed for future stud- 
ies. 

We finally comment that the case of dynam- 
ical, periodic, adjoint fermions was treated an- 
alytically in [12 , where it was shown that the 
chiral symmetry is indeed restored. 

7. Short-range correlations of dyons in 
quenched and non-quenched ensembles 

At high temperature, the vacuum is mostly 
dominated by perturbative fluctuations: the 
coupling is simply to small to allow any kind 
of large quantum effects. As we lower the tem- 
perature, the formation of topological objects 
starts to be possible. Although still suppressed, 
the vacuum is able to polarize into topological 
objects, which can support localized fermionic 
modes of small eigenvalue, but still not small 
enough to break chiral symmetry. Also, the 
"mass" M of these zero modes, which interpo- 
lates from 27rT to ttT i.e. reduces by a half 
in units of temperature as temperature is de- 
creased, makes it harder to break chiral sym- 
metry at high temperature, as the tail of zero 
modes, being the inverse of this mass, does not 
extend very far, and, thus, at small density of 
topological charge, the off-diagonal matrix el- 
ements are simply too small to matter, and 
the ensemble of dyons is the ensemble of LL 
pairs and neutral random clouds of light M, M 
dyons. 

However things change drastically as temper- 
ature is decreased. The gluon dynamics facili- 
tate the increase of topological density, due to 
a suppression of the action of any field config- 
uration at lower temperature by the coupling 
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1/^^. At one point the moduli space metric 
of topological objects becomes the sole dicta- 
tor of the distribution of topological (dyonic) 
field configuration of dyons, the rest simply be- 
ing fiuctuations which is also not suppressed at 
all, but does not change the important topo- 
logical properties of the background configura- 
tion. The dynamics of dyons becomes impor- 
tant solely through the geometry of the moduli 
space metric. 

The interaction of L and anti-L is very similar 
to the interaction between L and L, i.e. they 
repel (Coulombic-like), and the L and anti-M 
attract from the point of view of the metric (at 
large distances), just like L and M. The vacuum 
then needs to undergo a transition in the struc- 
ture, from pairs of LL and MM^ to the crystal 
of alternating dyons and antidyons. 

Such an abrupt change in the vacuum struc- 
ture of the quenched ensemble is absent upon 
introduction of dynamical fermions. In this 
case the situation changes drastically. The in- 
crease of topological density is suppressed by 
the presence of fermions, as increasing den- 
sity means making pairs come closer together. 
By arguments in the article, such a scenario 
will make the fermionic determinant be smaller 
and smaller, eventually going to zero if the 
molecules overlap. That means that a the same 
topological density is harder to develop with 
fermions. Increasing the fiavors makes it even 
harder, as the smallness of the fermionic mea- 
sure is enhanced by Nf. However, if Nf is not 
too large chiral symmetry may still be broken, 
but the nature of this transition is now vastly 
different. The pairs, instead of abruptly chang- 
ing their structure, are acting like slippery ob- 
jects, trying to keep their distance as large as 
possible to all other pairs. 



on the lattice, in which some of the tests 
proposed above were successfully performed. 
For periodic fermions their near-zero eigen- 
modes with small eigenvalues are indeed iden- 
tified with the type-M, M dyons, while for 
antiperiodic fermions those are indeed related 
with the type-L, L dyons. This conclusion was 
reached by correlating the topology with the 
sign of the Polyakov line at the dyon center. 
The difference in inverse participation ratios 
between the two cases is indeed naturally ex- 
plained by large difference in M and L ac- 
tions. For the "heavier" L dyons, even the 
shape of the fermionic eigenmodes was shown 
to agree with the corresponding semiclassical 
predictions. 
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Appendix A: Exact solution of the zero 
modes 



Although zero modes of a caloron were found 
in all generality elsewhere [H [10] , we here use 
an approach which is more illuminating. One of 
us would like to thank A. G. Abanov for useful 
discussions on this issue. 



Here we solve equations (32), where 



8. Small comment about deconfiement 

Recent work by Bruckmann et al [42 have ar- 
gued that uncorrelated dyons of all kinds create 
linear confinement. An ensemble of correlated 
neutral "molecules" we propose would not do 
that and generate short-range correlations only. 
This potentially links both deconfinement and 
chiral restoration with molecule formation. 

Note added: After our paper was com- 
pleted we learned about important study [49] 



r 

^ ^ l-vr/smh{vr) 
r ' 

We will take the lower sign (antiselfdual solu- 
tion). To do this we separate the matrix M(r) 
as 

M(r) = Mo(r) + Mi(r) , (A3) 
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where 



This equation has a general solution 



Mo{r) = (I + I) 1 



1 -\- rv coth(r'u) 



2r 



Mi(r) = 



^--1^3 
(^3 . 



sinh('L'r) 



(A4) 



(A5) 



The solution can then be written as a 
lo ^o{r)dr)x, or 



1 



a 



V r sinh 



rv 



(A6) 



Xi(0 = ci (^-2? + tanh e'^«+C2(2<r+tanh |)e 

(A13) 

with arbitrary constants ci,2- Using the first 
order equations we can write X2 as 

X2(0 = ci (^2? - coth e'^«+C2(2<r+coth |)e-^ 

(AM) 

The function X2(0 divergent when ^ ^ 0, ex- 
cept if ci = C2, in which case ^2(0) = 0- There- 
fore C2 = ci. The constant ci can be deter- 
mined by overall normalization. The solution 
then becomes 



Note that if we took the upper sign in eq. ( Al ), 



we would get a factor of ^/sinh(vr) in front of 
the solution. This is clearly non-normalizable, 
as it should be by the index theorem. 
The differential equation for x is 



-Mi(r)x, 



I.e. 



Xi{r) 
X2{r) 



smh('yr) p 
smh{vr) p 



(A7) 

(A8) 
(A9) 



we may take a change of variables ^ = rv. Then 
the equations read 



1 



:Xi{r) - TX2(r) , 



(AlO) 



sinh(^) ' 

X2(e) = -^^X2{r) - <;xi{r) , (All) 

where we labeled ^ = (j)/{yP) We now elimi- 
nate ^2, and obtain the second order differential 
equation 



2X1 



2 cosh^ I 



Xi 



-exi 



(A12) 



Xi(0 = 2ci ( -2(^sinh(^(^) + tanh - cosh(^(^; 



(A15a) 



X2(0 = 2ci ( 2c;cosh(^c;) - coth - sinh(^(^) 



(A15b) 



Finally, combining with ( A6 ) we obtain 



Xl,2 



(A16) 



^/v can be absorbed into constant ci, and our 
final expression is 



<^l,2 



Xl,2 



(A17) 



with functions xi,2 given by (|A6|), ^ = vr^ ^ = 
(j)/{yP). Note that the value of 0) is 

given by 



ci(l-4^') 



(A18) 



and the solution is completely regular at r = 0. 
Remarkably it turnes out that for ci = 1/2, the 
solution is already normalized (in the sense of 
J d^^^ integration). 



[1] F. Bruckmann, D. Nogradi and P. van arXiv:hep-th/0309008 . 

Baal, Acta Phys. Polon. B 34, 5717 (2003) [2] D. Diakonov, Topology and Confinement^ 



31 



arXiv:0906.2456V l (2009) 



E. M. Ilgenfritz and E. V. Shuryak, Nucl. Phys. 
B 319, 511 (1989). 

F. Bmckmann, PoS CONFINEMENTS, 179 
(2008) arXiv:090L0987 [hep-ph]]. 

J. Liao and E. Shuryak, Phys. Rev. C 75, 

054907 (2007) arXiv:hep-ph /061113l . 

A. A. Belavin, A. M. Polyakov, A. S. Schwartz 

and Yu. S. Tyupkin, Phys. Lett. B 59, 85 

(1975). 

T. Schafer and E. V. Shuryak, Rev. Mod. Phys. 

70, 323 (1998) arXiv:h ep-ph/9610451 . 

K. -M. Lee an d C. -h. Lu, Phys. Rev. D 58, 

025011 (1998) hep-th/9802108 . 

K. M. Lee and P. Yi, Phys. Rev. D 56, 3711 

(1997) arXiv:hep-th/9702107 . 

T. C. Kraan and P. van Baal, Phys. Lett. B 

435, 389 (1998) arXiv:hep-th/9806034 . 

M. Unsal, Phys. Rev. D 80, 065001 (2009) 

arXiv:0709.3269 [hep-th]]. 
M. Unsal, Phys. Rev. Lett. 100, 032005 (2008) 

arXiv:0708.1772 [hep-th]]. 
K. Miura, M. P. Lombardo and E. Pallante, 
arXiv:1110.3152 [hep-lat]. 

A. Deuzeman, M. P. Lombardo and E. Pal- 
lante, arXiv: 1012.6023 [hep-lat]. 
A. Cheng, A. Hasenfratz and D. Schaich, 
arXiv:1111.2317 [hep-lat]. 

G. Cossu, M. D'Elia, A. Di Giacomo, 
G. Lacagnina and C. Pica, Phys. Rev. D 77, 
074506 (2008) arXiv:0802.17 95 [hep-lat]]. 

R. D. Mawhinney, Nucl. Phys. Proc. Suppl. 
83, 57 (2000) arXiv:hep-lat/0001032 . 
X. Y. Jin and R. D. Mawhinney, PoS 
LAT2009, 049 (2009) arXiv:0910.3216| [hep- 
lat]]. 

F. Bruckmann, D. Nogradi and P. van Baal, 
_Nucl. Phys. B 666, 197 (2003) ar Xiv:hep- 
th/0305063. 

E. Bilgici, F. Bruckmann, J. Danzer, C. Gat- 
t ringer, C. Hagen, E. M. Ilgenfritz and 
A. Maas, Few Body Syst. 47, 125 (2009) 
arXiv:0906.3957 [hep-lat]]. 

E. Bilgici, C. Gattringer, E. -M. Ilgen- 
fritz and A. Maas, JHEP 0911, 035 (2009) 
arXiv:0904.3450 [hep-lat]]. 

C. Gattringer, M. Gockeler, P. E. L. Rakow, 
A. Schafer, W. Soldner and T. W ettig, Nuc l. 
Phys. Proc. S uppl. 106, 492 (2002) |arXiv:hep-| 
Jat/0110182^ 
C. Gattringer and S. Schaefer, Nucl. Phys. B 
654, 30 (2003) arXiv:hep-lat/0212029 . 
C. Gattringer, Phys. Rev. Lett. 88, 221601 
(2002) |arXiv:hep-l at/0202002 . 

F. Bruckmann, T. G. Kovacs and S. Schieren- 
berg, Phys. Rev. D 84, 034505 (2011) 



[26 

[27] 



arXiv: 1105.5336] [hep-lat]] . 
S. Hands and P. Kenny, Phys. Lett. B 701, 373 
(2011) arXiv:1104.0522 [hep-lat]]. 
M. N. Chernodub and V. I. Zakharov, Phys. 
Rev. L ett. 98, 082002 (2007) arXivihep-l 



ph/0611228 , 



[28] J. Liao and E. Shuryak, Phys. Rev. Lett. 101, 
162302 (2008) arXiv:0804.0255 [hep-ph]]. 

[29] H. Y. Chen, N. Dorey and K. Petunin, JHEP 
1006, 024 (2010) arXiv: 1004.0703 [hep-th]]. 

[30] R. Jackiw and C. Rebbi, Phys. Rev. D 13, 3398 
(1976). 

[31] Y. .M. Shnir, Berlin, Germany: Springer 
(2005) 532 p 

[32] P. van Baal, PoS C0NFINEMENT8, 062 

(2008) arXiv:0901.2853 [hep-ph]]. 

[33] E. V. Shuryak, Sov. Phys. JETP 47, 212 

(1978) [Zh. Eksp. Teor. Fiz. 74, 408 (1978)]. 
[34] R. D. Pisarski and L. G. Yaffe, Phys. Lett. B 

97, 110 (1980). 
[35] E. V. Shuryak, Nucl. Phys. B 203, 140 (1982). 
[36] D. Diakonov, N. Gromov, V. Petrov and 

S. Slizovskiy, Phys. Rev. D 70, 036003 (2004) 

arXiv:hep-th/0404042 . 
[37] E. Poppitz, M. Unsal, JHEP 1107, 082 (2011). 

arXiv:1105.3969 [hep-th]]. 
[38] E. V. Shuryak, arXiv:hep-ph/9909458| 
[39] E. Poppitz and M. Unsal, JHEP 0903, 027 

(2009) arXiv:0812.2085 [hep-th]]. 

[40] M. Garcia Perez, A. Gonzalez- Arroyo 
and A. Sastre, JHEP 0906, 065 (2009) 
[arXiv:0905.0645 [hep-th]]. 

[41] S. Borsanyi, Z. Fodor, C. Hoelbhng, 
S. D. Katz, S. Krieg, C. Ratti and K. K. Szabo 
[Wuppertal-Budapest Collaboration] , JHEP 
1009, 073 (2010) a rXiv: 1005.3 508 [hep-lat]]. 

[42] F. Bruckmann, S. Dinter, E. M. Ilgenfritz, 
B. Maier, M. MuUer-Preussker and M. Wag- 
ner, arXiv: 111 1.3158 [hep-ph]. 

[43] K. -Y. Kim, S. -J. Sin and I. Zahed, JHEP 
0809 (2008) 001 arXiv:0712.1582 [hep-th]]. 

[44] C. Ratti and E. Shuryak, Phys. Rev. D 80, 
034004 (2009) arXiv:0811.4174 [hep-ph]]. 

[45] N. S. Manton, Nucl. Phys. B 126, 525 (1977). 

[46] J. N. Goldberg, P. S. Jang, S. Y. Park and 
K. C. Wah, Phys. Rev. D 18, 542 (1978). 

[47] S. F. Magruder, Phys. Rev. D 17, 3257 (1978). 

[48] S. Gabor. Orthogonal Polynomials, American 
Mathematical Society 1995. 

[49] V. G. Bornyakov, E. -M. Ilgenfritz, 
B. V. Martemyanov and M. MuUer- 
Preussker, Phys. Rev. D 79, 034506 (2009) 
arXiv:0809.2142 [hep-lat]]. 

[50] In the text we often refer to the "holonomy 
mass" or simply "mass", occasionally with- 
out the quotation marks, to indicate the mass 



32 



given by the Higgs effect of A4. This "mass" 
does not break chiral symmetry and should in 
no way be confused with the fermion mass. In 
our study fermions are taken as massless. 
[51] for expUcit form of function 7i,A see the Ap- 
pendix 

[52] This ansatz is inspired by the fact that the 
fermions have two SU{2) indices, one spin and 
one color. Since we consider the fundamen- 
tal color representation (i.e. both are 2 rep- 
resentations), these indices can combine into 
2 2 = 3, which is precisely the decompo- 
sition considered in the text. See e.g. 31 

[53] In 1977, when QCD was only 4 years old, 
the main finding was its positivity, which en- 



sured screening of a charge, as opposed to the 
antiscreening by the vacuum loops: thus the 
"plasma" name. 

[54] Note that he explicit numerical factors and the 
definition of the A parameter are not really 
important since they will change in a theory 
with fermions. 

[55] We also ignored small nondiagonal matrix ele- 
ments of the type m < Ii\l2 > resulting from 
non-orthogonality of the zero modes of differ- 
ent pseudoparticles. 

[56] This dip should not be confused with that due 
to the finite size of the system, which is well 
described by the random matrix theory. 



